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Theoretical  methods  are  reported  for  ab  initio  calculations  of  the  adiabatic  (Born-Oppenheimer) 
electronic  wave  functions  and  potential  energy  surfaces  of  molecules  and  other  atomic  aggregates. 

An  outer  product  of  complete  sets  of  atomic  eigenstates  familiar  from  perturbation-theoretical 
treatments  of  long-range  interactions  is  employed  as  a  representational  basis  without  prior 
enforcement  of  aggregate  wave  function  antisymmetry.  The  nature  and  attributes  of  this  atomic 
spectral-product  basis  are  indicated,  completeness  proofs  for  representation  of  antisymmetric  states 
provided,  convergence  of  Schrodinger  eigenstates  in  the  basis  established,  and  strategies  for 
computational  implemention  of  the  theory  described.  A  diabaticlike  Hamiltonian  matrix 
representative  is  obtained,  which  is  additive  in  atomic-energy  and  pairwise-atomic 
interaction-energy  matrices,  providing  a  basis  for  molecular  calculations  in  terms  of  the 
(Coulombic)  interactions  of  the  atomic  constituents.  The  spectral-product  basis  is  shown  to  contain 
the  totally  antisymmetric  irreducible  representation  of  the  symmetric  group  of  aggregate  electron 
coordinate  permutations  once  and  only  once,  but  to  also  span  other  (non-Pauli)  symmetric  group 
representations  known  to  contain  unphysical  discrete  states  and  associated  continua  in  which  the 
physically  significant  Schrodinger  eigenstates  are  generally  embedded.  These  unphysical 
representations  are  avoided  by  isolating  the  physical  block  of  the  Hamiltonian  matrix  with  a  unitary 
transformation  obtained  from  the  metric  matrix  of  the  explicitly  antisymmetrized  spectral-product 
basis.  A  formal  proof  of  convergence  is  given  in  the  limit  of  spectral  closure  to  wave  functions  and 
energy  surfaces  obtained  employing  conventional  prior  antisymmetrization,  but  determined  without 
repeated  calculations  of  Hamiltonian  matrix  elements  as  integrals  over  explicitly  antisymmetric 
aggregate  basis  states.  Computational  implementations  of  the  theory  employ  efficient  recursive 
methods  which  avoid  explicit  construction  the  metric  matrix  and  do  not  require  storage  of  the  full 
Hamiltonian  matrix  to  isolate  the  antisymmetric  subspace  of  the  spectral-product  representation. 

Calculations  of  the  lowest-lying  singlet  and  triplet  electronic  states  of  the  covalent  electron  pair 
bond  (H2)  illustrate  the  various  theorems  devised  and  demonstrate  the  degree  of  convergence 
achieved  to  values  obtained  employing  conventional  prior  antisymmetrization.  Concluding  remarks 
place  the  atomic  spectral-product  development  in  the  context  of  currently  employed  approaches  for 
ab  initio  construction  of  adiabatic  electronic  eigenfunctions  and  potential  energy  surfaces,  provide 
comparisons  with  earlier  related  approaches,  and  indicate  prospects  for  more  general  applications  of 
the  method.  ©  2004  American  Institute  of  Physics.  [DOI:  10.1063/1.1794634] 


I.  INTRODUCTION 

Adiabatic  electronic  wave  functions  provide  well-known 
first  approximations  to  the  physical  attributes  of  molecules 
and  other  forms  of  matter, 1-3  and  are  prerequisites  for  studies 
of  the  nonadiabatic  dynamical  consequences  of  avoided 
crossings  and  conical  intersections  in  the  potential  energy 
surfaces  that  guide  the  pathways  of  chemical  reactions.4-6 
Methods  currently  employed  for  constructing  such  wave 
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functions  and  energy  surfaces  include  ab  initio  quantum 
chemistry,7-10  quantum  Monte-Carlo  simulations,11,12 
density-functional  approaches,13-18  semiempirical  techni¬ 
ques,19-2'  hybrid  quantum  mechanical/molecular  mechanical 
(QM/MM)  combinations,24  and  various  forms  of  long-range 
perturbation  theory.25-38  Although  these  methods  have  pro¬ 
vided  a  significant  amount  of  useful  information  over  the 
past  few  decades,  they  arguably  do  not  constitute  a  quantum 
theory  of  chemical  bonding,  which  would  preferably  express 
the  chemical  attributes  of  compounds  in  terms  of  intrinsic 
properties  of  their  atomic  constitutents.  New  approaches  in 
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this  spirit  to  construction  of  adiabatic  electronic  wave  func¬ 
tions  are  clearly  welcome,  particularly  if  they  can  retain  the 
desirable  features  of  the  existing  methods  while  ameliorating 
some  of  their  shortcomings. 

The  ab  initio  approaches  employing  multiconfigurational 
Hartree-Fock,  Mpller-Plessett  perturbation  theory,  and  mul¬ 
tireference  configuration  mixing  methods,7'8  and  the 
coupled-cluster,9  valence-bond,10  and  quantum  Monte 
Carlo1112  approaches,  can  provide  highly  accurate  energies 
and  other  properties  for  the  ground  and  excited  states  of 
small  molecules,  but  they  generally  require  allocations  of 
significant  computational  resources  for  this  purpose, 
and  they  are  apparently  not  yet  applicable  to  large  atomic 
aggregates.  By  contrast,  the  widely  employed 
density-functional, 13-18  semiempirical,19  21  and  QM/MM24 
approaches  can  be  applied  to  larger  aggregates,  including 
atomic  clusters,22  condensed  matter  systems,23  and  biological 
macromolecules,24  and  in  certain  cases  provide  electronically 
excited  states,  although  their  a  priori  accuracy  is  frequently 
unknown.  Long-range  Rayleigh-Schrodinger  and  symmetry- 
adapted  perturbation  theories  have  also  proved  useful  in  ap¬ 
propriate  limits,25  38  although  widespread  and  systematic  ap¬ 
plications  of  perturbation  theory  to  chemical  interactions  are 
apparently  not  yet  in  evidence.  These  methods  have  all  con¬ 
tributed  to  the  collective  development  of  a  significant  metho- 
logical  and  computational  base  from  which  to  draw  upon  in 
devising  new  ab  initio  approaches  to  electronic-structure  and 
spectral  calculations. 39-44 

In  the  present  paper,  theoretical  methods  are  described 
for  ab  initio  calculations  of  the  adiabatic  electronic  wave 
functions  and  potential  energy  surfaces  of  molecules  and 
other  atomic  aggregates.45  Certain  of  the  desirable  attributes 
of  the  computational  methods  currently  in  play  for  these  pur¬ 
poses  are  incorporated  in  the  development,  which  employs 
an  outer  product  of  complete  sets  of  atomic  spectral  eigen¬ 
functions  familiar  from  the  aforementioned  perturbation- 
theoretical  treatments  of  long-range  atomic  interactions.25  In 
contrast  to  related  earlier  developments  which  employ  prior 
antisymmetrization  of  products  of  atomic  wave  functions  in 
atoms-in-molecules  approaches,46-58  antisymmetry  is  en¬ 
forced  in  the  present  development  subsequent  to  construction 
of  the  Hamiltonian  matrix  employing  a  unitary  transforma¬ 
tion  to  isolate  its  antisymmetric  subspace.  The  total  elec¬ 
tronic  Hamiltonian  operator  in  the  atomic-product  represen¬ 
tation  takes  a  matrix  form  in  which  the  individual  atomic 
energies  appear  explicitly  as  constants  independent  of  atomic 
positions,  with  only  a  sum  of  pairwise  atomic  interaction- 
energy  matrices  requiring  repeated  evaluation  as  a  function 
of  aggregate  atomic  configuration  in  potential-energy-surface 
calculations.  The  pairwise  interaction-energy  matrices  de¬ 
pend  only  upon  the  vector  separations  of  the  two  interacting 
atoms  and  are  independent  of  the  positions  of  all  other  atoms 
in  the  aggregate,  and  can  be  expressed  in  terms  of  one- 
electron  atomic  charge  and  transition  densities,  which  pro¬ 
vide  the  computational  invariants  of  the  development.  Ac¬ 
cordingly,  the  present  approach  does  not  require  repeated 
calculations  of  Hamiltonian  matrix  elements  as  integrals  over 
antisymmetric  many-electron  aggregate  basis  functions; 
rather,  one-electron  atomic  charge  and  transition  densities  are 


calculated  once  and  for  all  for  a  spectrum  of  atomic  states 
and  retained  for  repeated  use  in  construction  of  the  Hamil¬ 
tonian  matrix  for  arbitrary  spatial  arrangements  of  interacting 
atoms. 

Notational  conventions  are  established  and  the  (Coulom- 
bic)  form  of  the  many-electron  Hamiltonian  operator  em¬ 
ployed  in  the  development  given  in  Sec.  II.  The  atomic 
spectral-product  basis  and  certain  of  its  attributes  are  de¬ 
scribed  in  Sec.  Ill  employing  complete  sets  of  spin  orbitals  in 
representations  of  the  many-electron  atomic  spectral  eigen¬ 
states.  It  is  shown  that  the  totally  antisymmetric  irreducible 
representation  of  the  aggregate  symmetric  group  is  spanned 
once  and  only  once  by  the  atomic -product  basis,  although 
other  (non-Pauli)  symmetric -group  representations  are  also 
contained  in  the  reducible  product.  The  form  of  the  Hamil¬ 
tonian  matrix  in  the  spectral-product  representation  is  given 
in  Sec.  IV,  methods  are  described  for  its  evaluation  in  terms 
of  one-electron  atomic  transition  density  matrices,  and  cer¬ 
tain  attributes  of  the  spectrum  of  the  Hamiltonian  matrix  are 
indicated.  In  Sec.  V,  methods  are  reported  for  formal  isola¬ 
tion  of  the  totally  antisymmetric  block  of  the  spectral- 
product  Hamiltonian  matrix  employing  a  unitary  transforma¬ 
tion  derived  from  the  metric  matrix  of  the  antisymmetrized 
spectral-product  basis.  Convergence  is  demonstrated  in  this 
way  to  eigenstates  and  energies  identical  to  those  obtained  in 
the  linearly  independent  subspace  of  this  prior  antisymme¬ 
trized  representation.  Aspects  of  the  theory  are  illustrated 
with  calculations  reported  in  Sec.  VI,  which  employ  efficient 
recursive  methods  for  isolating  the  totally  antisymmetric 
block  of  the  Hamiltonian  matrix,  avoiding  explicit  construc¬ 
tion  of  the  metric  matrix  and  of  the  entire  Hamiltonian  ma¬ 
trix  at  any  one  time.  Convergence  to  the  lowest-lying  singlet 
and  triplet  potential  energy  curves  of  the  electron  pair  bond 
(H2)  is  achieved  employing  even-tempered  Gaussian  orbital 
representations  and  related  computational  methodology.  Fi¬ 
nally,  the  development  is  discussed  in  the  context  of  cur¬ 
rently  employed  electronic-structure  approaches,  compari¬ 
sons  are  made  with  related  earlier  theoretical  methods,  and 
prospects  for  more  general  applications  of  the  method  are 
indicated  in  concluding  remarks  provided  in  Sec.  VII. 

II.  THE  AGGREGATE  SCHRODINGER  EQUATION 

The  adiabatic  electronic  eigenstates  and  energies  of  an 
/V-atom  molecule  or  other  atomic  aggregate  are  obtained 
from  the  Schrodinger  equation 

//(l,2,...,n:R)'I'(l,2,...,n:R)  =  'E(l,2,...,n:R)-E(R),  (1) 

where  M,(l,2,...,n:R)  is  a  row  vector  of  orthonormal  eigen¬ 
states  AFr(L2,...,n:R)  having  quantum  labels  T  and  E(R)  is 
the  diagonal  matrix  of  corresponding  eigenvalues  £r(R). 
The  Hamiltonian  operator  appearing  in  Eq.  (1)  is  written  in 
the  form 

//(1,2,  ...,n:R) 

N  (  N  I 

=  2  U(a)(i:RJ+  2  V(“’«(i;j:R^)  ,  (2) 

a=  1  y  /3=  a+  1 

where 
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H{a)(  i:Ra)  =  2 
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is  the  Hamiltonian  operator  in  the  familar  Coulomb  approxi¬ 
mation  for  a  set  of  n  a  electrons  which  are  arbitrarily  associ¬ 
ated  with  the  atomic  nucleus  at  position  R„ ,  and  whose  3  na 
spatial  coordinates  are  represented  by  the  single  vector  label 
i,  and 


(i;j:Raie) 


7  rr  2  a  7  2  ^  R  r~7  2 

Ra/3  i  rif3  j  rja 


"a  n(3  .2 

+2  2—  (4) 

<■  j 


is  the  Coulombic  interaction  potential  between  atoms  a  and 
/3,  where  the  coordinate  label  j  refers  to  the  3  rip  spatial  elec¬ 
tron  coordinates  arbitrarily  associated  with  the  atomic 
nucleus  at  the  position  Rg.59  The  indecies  1,  2,...,n  in  Eqs. 
(1)  and  (2)  represent  the  collected  spatial  coordinates  of  dis¬ 
joint  sets  of  electrons,  each  arbitrarily  associated  with  one  of 
the  atomic  nuclei  in  the  aggregate. 

In  the  foregoing  equations  and  throughout  the  develop¬ 
ment,  all  electron  (r;)  and  atomic -position  (Ra)  coordinates 
are  defined  in  the  laboratory  frame,  r/a=r,  — Ra  and  r//3 
—  Yj —  Rg  are  electron  coordinate  vectors  relative  to  the  indi¬ 
cated  atomic  positions,  r(/  =  r(  — r;  is  the  (i,j)  electron  vector 
separation,  Rajg=Ra—  R^  is  the  (a,/3)  atomic-position  vector 
separation,  R=(R1,R2,...,R;v)  specifies  the  positions  of  all 
the  atoms  in  the  entire  aggregate,  e  is  the  magnitude  of  the 
electronic  charge,  m  is  the  electron  mass,  Za  and  Zp  are 
atomic  numbers,  and  the  other  symbols  have  their  usual 
meanings.59  A  colon  is  employed  to  separate  electron  coor¬ 
dinates  from  atomic  positions  (i:Ra),  and  a  semicolon  is 
employed  in  Eq.  (4)  to  distinguish  sets  of  electrons  (i;j)  be¬ 
tween  which  the  interaction  potential  V/(“+^(i;j:  Rap)  is  not 
totally  symmetric  under  coordinate  permutations.  Of  course, 
the  n,  electron  (nt=nl+n2^ - aggregate  Hamil¬ 

tonian  //(l,2, ...,n:R)  and  the  individual  n ^-electron  atomic 
Hamiltonians  H(a)(i:R  a )  are  symmetric  sums  of  one-  and 
two-electron  operators  in  their  respective  electron  labels.  Ac¬ 
cordingly,  the  notation  employed  in  Eqs.  (1)— (4)  does  not  by 
itself  constitute  a  meaningful  assignment  of  particular  elec¬ 
trons  to  particular  atoms,  although  it  does  suggest  this  pos¬ 
sibility  in  conjunction  with  the  choice  of  an  appropriate  elec¬ 
tronic  representational  basis. 

The  physical  eigenstates  'P'<p'l(l,2,...,  n:R)  and  eigenval¬ 
ues  E(p\R)  included  in  the  solutions  of  Eq.  (1)  for  the 
Hamiltonian  of  Eqs.  (2) — (4)  have  good  total  spin  and  spin- 
projection  quantum  numbers  ( S,Ms )  and  transform  as  the 
totally  antisymmetric  representation  (1"')  of  the  symmetric 
group  S n<  of  the  aggregate  electron  spin  and  spatial  coordi¬ 
nate  permutations,60  whereas  additional  (non-Pauli)  eigenso- 
lutions  *P(1|)(1,2,. .  .,n:R)  of  Eq.  (1)  which  transform  as  other 
than  the  totally  antisymmetric  irreducible  representation  un¬ 
der  electron  coordinate  permutations  are  generally  regarded 
as  unphysical,  as  are  the  corresponding  energies  E(i,)(R).61 


III.  DEFINITION  AND  ATTRIBUTES 
OF  THE  SPECTRAL-PRODUCT  BASIS 

The  physical  and  unphysical  solutions  of  Eq.  (1)  which 
together  comprise  the  full  spectrum  of  the  Hamiltonian  op¬ 
erator  of  Eqs.  (2)  to  (4)  can  be  constructed  in  a  many- 
electron  spectral-product  basis  the  nature  and  attributes  of 
which  are  described  in  this  section.  The  presence  of  unphysi¬ 
cal  states  and  the  importance  of  devising  methods  for  their 
elimination  is  emphasized  throughout  the  development. 

A.  Spin-orbital  basis  sets 

Conventional  variational  approximations  to  the  solutions 
of  Eq.  (1)  in  their  most  common  form  employ  one-electron 
basis  sets  of  atomic  spin  orbitals  in  Slater,  Gaussian,  or  other 
forms  centered  at  each  of  the  aggregate  atoms  (a 
=  1,2 Such  basis  sets,  employed  here  largely  for  pur¬ 
poses  of  analysis,62  are  conveniently  written  as  row  vectors 
<f>{aXi\ R„)  of  orthonormal  one-electron  functions 
4la\r-Ra)  of  the  individual  electronic  spin  and  of  spatial 

'a 

coordinates  measured  relative  to  the  atomic  positions  ( i  =  r, 
—  Ra).  These  functions  generally  enumerated  employing  an 
appropriate  set  of  quantum  labels  ya  are  assumed  complete 
in  a  finite  three-dimensional  spatial  domain,  and,  accord¬ 
ingly,  satisfy  the  closure  relation 

00 

^a\v.Ra)-^a\i':Ra)^  2  Ra)4>\a\i'-Ra)* 

V  =1  r<*  ra 

/ a 

(5) 

in  this  domain,  where  </){a\i':Ra) 1  is  a  column  vector 
comprised  of  the  complex  conjugate  spin  orbitals 
<//“  (; ':RA*.  The  choice  of  particular  spin  orbitals  em- 

'a 

ployed  in  Eq.  (5)  is  arbitrary  in  the  limit  of  a  complete  set, 
but  the  rate  of  convergence  of  computational  developments 
built  upon  them  will  generally  depend  upon  the  specific  basis 
employed.62 

Traditional  uses  of  these  one-electron  basis  sets  in  con¬ 
nection  with  solution  of  Eq.  (1)  follow  well-known 
molecular-orbital  or  valence-bond  methods  in  constructing 
many-electron  representational  basis  sets.63  These  conven¬ 
tional  developments  have  in  common  the  representation  of 
each  electron  in  a  “symmetrical”  or  “democratic”  manner 
and  the  use  of  prior  term-by-term  electron  coordinate  anti- 
symmetrization  of  the  many-electron  basis  employed.61  In 
the  molecular-orbital  approaches,  the  common  one-electron 
basis  set  that  extends  over  the  3  n  ,-dimensional  spatial  do¬ 
main  of  the  aggregate  is  employed  to  describe  the  positions 
of  any  of  the  electrons  in  antisymmetric  Slater 
determinants,7,8  whereas  in  the  valence-bond  methods  over¬ 
all  term-by-term  prior  antisymmetry  ensures  the  canonical 
spin-paired  product  states  employed  treat  the  electrons  in  a 
democratic  fashion.10  Of  course,  use  of  prior  antisymmetry, 
and  of  symmetry  adaptation  more  generally,  guarantees  that 
the  matrix  representative  of  the  Hamiltonian  operator  ob¬ 
tained  will  be  block  diagonal  and  of  minimal  possible  dimen¬ 
sion  for  the  states  of  interest  in  the  basis  employed,  in  accor¬ 
dance  with  Schur’s  lemmas.60 


Downloaded  01  Nov  2007  to  198.14.54.253.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://jcp.aip.org/jcp/copyright.jsp 


9326 


J.  Chem.  Phys.,  Vol.  121,  No.  19,  15  November  2004 


Langhoff  et  al. 


B.  Many-electron  outer-product  basis  set 

In  contrast  to  the  aforementioned  conventional  develop¬ 
ments,  it  is  possible  to  forego  prior  symmetry  adaptation  of 
the  individual  terms  in  the  many-electron  basis  set  in  favor 
of  other  potential  advantages  and  simplifications  by  employ¬ 
ing  n  ,-electron  states  which  are  not  individually  antisymmet¬ 
ric  at  the  outset.  Such  states  can  be  constructed  from  prod¬ 
ucts  of  n  ^-electron  atomic  basis  states  which  themselves  are 
written  in  the  form  of  simple  products  of  spin  orbitals, 

nc 

n^(i:Ra)  =  n  4%\ i:RJ,  (6) 

a  y  a 

"a 

where  the  label  specifies  a  particular  choice  of  spin  or¬ 
bitals  on  the  atom  a  in  the  indicated  n  ,rterm  product.  The 
row  vector  I  l(")(i:  Ra)  obtained  from  all  Ta  values  includes 
all  possible  unrestricted  spin-orbital  products  in  Eq.  (6), 
thereby  providing  a  complete  representational  basis  for 
n  ^-electron  eigenstates  consequent  of  Eq.  (5).  Accordingly, 
the  complete  outer  product  of  such  states  for  all  atoms  in  the 
aggregate,  written  in  the  form  of  the  complete  separable  Hil¬ 
bert  space, 

II(l;2;...;n:R) 

=  {n(1)(l:Ri)®n(2)(2:R2)®---®n(A')(n:RA,)}0,  (7) 

similarly  provides  a  complete  representational  basis  for  the 
entire  n  ,-electron  aggregate.  Here,  semicolons  are  employed 
to  separate  groups  of  electrons  (l;2;...;n)  distinguished  by 
virtue  of  their  representations  in  the  different  sets  of  spin- 
orbital  product  states  employed  for  each  atom,  whereas  the 
outer-product  symbol  “®”  is  used  to  indicate  that  all  pos¬ 
sible  products  of  the  atomic  basis  functions  are  included  in 
the  vector  II(l;2;...;n:R).  The  bracket  symbol  implies 

choice  of  a  particular  rule  for  enumerating  and  ordering  the 
product  functions  in  the  row  vector  II(l;2;...;n:R),  which 
rule,  although  arbitrary,  must  be  adhered  to  in  matrices  con¬ 
structed  in  the  outer-product  basis  and  in  related  vector  and 
matrix  multiplications. 

The  basis  of  Eq.  (7),  in  which  specific  sets  of  electrons 
(l;2;...;n)  are  arbitrarily  assigned  to  specific  atoms  (1,2 
and  are  described  employing  different  sets  of  orbital  func¬ 
tions  [Eq.  (5)],  does  not  provide  individual  n , -electron  func¬ 
tions  which  transform  irreducibly  under  aggregate  electron 
coordinate  permutations,  nor  are  these  functions  generally 
eigenstates  of  total  electron  spin  or  of  other  aggregate  com¬ 
muting  operators.  Nevertheless,  as  a  consequence  of  the  clo¬ 
sure  relations  of  Eq.  (5)  for  the  different  sets  of  spin  orbitals, 
the  aggregate  outer-product  basis  satisfies  the  closure  rela¬ 
tion 

n(l;2;...;n:R)  •  II(l';2';...;n':R)t 

N  N 

=  n  n(af)(i:Rj •  n(ff)(i' :Rjf^  II  <S(a)(i-i'),  (8) 

a= 1  a= 1 

where  the  product  of  S  functions  on  the  right-hand  side  ex¬ 
presses  completeness  of  the  basis  as  a  separable  Hilbert 
space  in  the  spin  coordinates  and  over  a  3  n, -dimensional 
spatial  domain.  Accordingly,  convergence  at  least  in  the 
norm  to  physical  and  unphysical  eigenstates  of  irreducible 


symmetry  can  be  achieved  in  the  basis  of  Eqs.  (5)-(8),  po¬ 
tentially  providing  approximations  to  any  of  the  aggregate 
Schrodinger  eigensolutions  of  Eq.  (1). 


C.  Atomic  spectral  eigenstates 

It  is  convenient  to  perform  a  series  of  formal  unitary 
transformations  of  the  complete  product  basis  sets  of  Eq.  (6) 
which  bring  the  matrix  representatives  of  the  individual 
atomic  Hamiltonian  operators  of  Eq.  (3)  into  diagonal  forms, 
thereby  simplifying  the  total  Hamiltonian  matrix  but  leaving 
the  solutions  of  Eq.  (1)  in  the  basis  invariant.  Specifically, 
because  a  single  spin-orbital  basis  set  is  used  for  all  electrons 
i  on  atom  a,  and  all  spin-orbital  products  are  present  in  Eq. 
(6),  all  permutations  among  electron  coordinates  i  are  also 
included  explicitly  for  each  atomic  site.  It  should  be  recalled 
in  this  connection  that  all  permutations  of  electron  coordi¬ 
nates  among  the  spin-orbital  products  of  Eq.  (6)  are  included 
explicitly  in  conventional  atomic  structure  theory,  so  that  the 
orbital  configurations  employed  to  represent  atomic  states 
must  be  limited  to  sets  of  ordered  spin-orbital  products  in 
Slater  determinants  to  insure  the  linear  independence  of  the 
basis  states  employed.62  The  present  development,  in  which 
all  orderings  of  spin-orbital  products  for  a  given  atom  are 
included  but  only  a  single  assignment  of  electrons  to  spin 
orbitals  is  made,  is  entirely  equivalent  to  the  standard 
method  because  a  common  spin-orbital  basis  is  employed  in 
a  symmetrical  representation  of  all  the  electrons  in  a  given 
atom.61 

The  eigenstates  of  the  atomic  Hamiltonian  of  Eq.  (3), 
H(a\i:Ra)&a\l.Ra)  =  &a\i:Ra)-Eia\  (9) 


are  represented  in  the  form 

<K(“)(i:RJ  =  n(“)(i:RJ-U(“),  (10) 

where  the  unitary  transformation  matrix  U(  n)  is  independent 
of  Ra,  and  E(a)  is  the  diagonal  matrix  of  all  atomic  eigen¬ 
values  of  the  atom  a  obtained  from  the  matrix  Schrodinger 
equation 


(n(“)(i:Ra)|H(  “)(i:Ra)|n<“)(i:RJ)-U(“) 


(ID 


In  Eq.  (11),  an  outer- vector-product  convention  is  employed 
on  the  left-hand  side  to  provide  matrices  having  conventional 
row  and  column  labels  formed  from  the  indicated  bra 
(n(“)(i:R„)|  and  ket  I ll'*',(i: R„))  vectors.  The  irreducible- 
symmetry  labels  or  good  quantum  numbers  of  the  eigenstates 
of  Eq.  (9)  can  also  be  accommodated  into  the  transformation 
matrix  of  Eqs.  (10)  and  (11)  without  loss  of  generality,  and 
so  0<a)(i:RJ  can  be  taken  to  be  a  row  vector  of  functions 
d4a)(i:RJ  which  are  multiplet  eigenstates  of  the  Hamil- 

1  a 

tonian  H{a\i:Ra)  having  specific  energy  Ey  ,  orbital  and 
spin  angular  momentum,  and  parity  values;  i.e.,  T„ 
=  (E,L,M L  ,S ,M s  ,P) a  ,63  These  atomic  eigenfunctions 
transform  either  as  physical  (totally  antisymmetric)  or  un¬ 
physical  (non-Pauli)  irreducible  representations  under  the  co¬ 
ordinate  permutation  group  Sn  for  electrons  in  the  set  i.60 
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Combining  the  unitary  transformations  of  Eqs.  (10)  and 
(11)  for  all  atoms  in  the  aggregate,  a  complete  set  of  products 
of  atomic  spectral  eigenstates  is  obtained  in  the  form 

<J>(l;2;...;n:R)  =  II(l;2;...;n:R)-U,  (12) 

where  U={U(1)®U(2)<8>- ■  -<g>Uw}0  is  the  appropriately  or¬ 
dered  outer  matrix  product  of  unitary  transformation  matri¬ 
ces  of  the  individual  spectral  basis  sets  of  Eqs.  (10)  and  (11), 
and 

<&(l;2;...;n:R) 

=  {<I»(1)(l:Ri)®<I»<2)(2:R2)®---®<I»(A')(n:Riv)}0  (13) 

is  the  resulting  complete  set  of  ordered  products  of  atomic 
spectral  eigenstates.  As  a  consequence  of  Eq.  (8)  and  the 
definitions  of  Eqs.  (9)— (12),  the  spectral-product  states  of 
Eq.  (13)  satisfy  the  closure  relation 

<I»(l;2;...;n:R)-«l>(l';2';...;n':R)t 

N  N 

=  11  0(“)(i:Rj.^“)(i':Ra)t^n  <S(a)(i-i')  (14) 

a= 1  a= 1 

for  all  irreducible  representations  of  Sn<,  and  so  can  be  em¬ 
ployed  in  descriptions  of  any  of  the  eigensolutions  of  Eq.  (1). 
Use  of  the  spectral -product  representation  of  Eq.  (13)  in 
place  of  the  outer-product  representation  of  Eq.  (7)  consti¬ 
tutes  a  partial  prediagonalization  of  the  aggregate  Hamil¬ 
tonian  matrix  which  leaves  its  dimension  unchanged  and  oth¬ 
erwise  provides  results  equivalent  to  those  obtained  from  the 
outer- product  basis. 

D.  Restriction  to  antisymmetric  atomic  states 

In  view  of  the  nonsymmetrical  nature  of  the  spectral- 
product  basis  of  Eq.  (13),  and  of  the  fact  that  the  outer  prod¬ 
uct  of  the  symmetric  groups  of  the  individual  atomic  con¬ 
stituents  forms  a  subgroup  of  the  full  electron-coordinate 
permutation  group  (Sl®S2®mmm®SNClSni),  some  care  must 
be  exercised  in  employing  standard  group-theoretical  proce¬ 
dures  in  identifying  and  isolating  the  desired  totally  antisym¬ 
metric  representation  (1"')  of  S„  contained  in  the  basis.  The 
appropriate  group-theoretical  product  in  this  case  is  the  outer 
product  of  subgroup  representations,  as  opposed  to  the  inner 
Kronecker  group  product  which  is  commonly  reduced  em¬ 
ploying  a  Clebsch-Gordan  decomposition.60  Moreover,  in 
the  representation  of  Eq.  (13),  although  electron-coordinate 
permutations  within  each  set  i  are  included  explicitly,  those 
permutations  involving  at  least  one  transposition  of  electron 
coordinates  from  different  sets  i  and  j  at  the  two  atomic  sites 
a  and  /3  are  spanned  implicitly.  Nevertheless,  in  the  limit  of 
the  closure  associated  with  use  of  all  orderings  of  spin- 
orbital  products  in  Eq.  (13),  the  latter  implicit  electron- 
coordinate  permutations  are  all  effectively  present  in  Eq. 
(13)  and  the  standard  outer-product  reduction  theory 
applies.60 

Although  there  is  apparently  no  convenient  way  in 
which  to  use  group-theoretical  methods  in  the  explicit  con¬ 
struction  of  a  particular  irreducible  subspace  in  the  ordered- 
product  representation  of  Eq.  (13),  the  outer-product  reduc¬ 
tion  procedure  can,  nevertheless,  be  employed  as  an 


analytical  device  in  this  connection.32  In  particular,  it  can  be 
demonstrated  that  some  (but  not  all)  of  the  unwanted  un¬ 
physical  representations  of  S„t  spanned  by  the  basis  of  Eq. 
(13)  are  removed  by  the  simple  expedient  of  restricting  the 
atomic  eigenstates  in  the  vectors  i:  Ra)  for  all  atoms  a 
to  the  physical  subspace  “p”  of  eigenstates  <1**,  (i:  Rff) 
which  are  antisymmetric  in  the  n  a  electron  coordinates  i.  The 
required  demonstration  follows  from  the  standard  outer- 
product  reduction  rules,60  which  indicate  that  the  presence  of 
any  horizontal  (“symmetric”)  strip  in  either  of  two  Young 
patterns  associated  with  distinct  symmetric  groups  disallows 
the  appearance  of  the  totally  antisymmetric  representation  in 
the  reduction  of  the  outer-product  of  the  two  patterns  in  the 
combined  permutation  group.  More  specifically,  the  reduc¬ 
tion  of  the  outer  product  of  any  two  irreducible  representa¬ 
tions  ( y„  <8>  ym)  of  two  distinct  symmetric  groups  S n ,  Sm  into 
irreducible  representations  of  the  group  Sn+m  includes  the 
totally  antisymmetric  representation  U'!+m>  of  Sn+m  once 
and  only  once  if  and  only  if  the  two  subgroup  representations 
are  the  totally  antisymmetric  representations  yn  =  1 "  and 
ym=  lm.  Letting  Sn  represent  the  symmetric  group  for  the 

first  a—  1  atoms  in  the  aggregate  (n  =  n\  +  n2-\ - 

setting  Sm  =  Sn  ,  the  symmetric  group  of  the  atom  a,  and 
assuming  that  1 "  is  contained  only  once  in  Sn ,  it  is  seen  that 
l"<8>y„  contains  1  (n+na>  oniy  once  if  and  only  if  y„ 
=  1"“.  Correspondingly,  it  follows  by  induction  that  the  to¬ 
tally  antisymmetric  representation  1"'  of  S„t  is  obtained  once 
and  only  once  in  the  reduction  of  the  entire  outer  product 
(y«j® y«2®‘ ■ -® ynN)  if  and  oniy  if  y«  =i"a  for  ail  a. 

In  accordance  with  the  foregoing  observations,  a  com¬ 
plete  but  not  overcomplete  linearly  independent  orthonormal 
subspace  of  the  spectral-product  basis  of  Eq.  (13)  suitable  for 
representations  of  totally  antisymmetric  aggregate  states  is 
contained  in  the  products  of  the  physical  (antisymmetric) 
atomic  eigenstates  ^“'(iiR^), 

<l>(l;2;...;n:R) 

=  {1>J,1)(l:R1)®^2)(2:R2)®-"®Of)(n:Rw)}0.  (15) 

The  spectral  closure  relation  of  Eq.  (14)  correspondingly 
now  becomes  operative  only  in  the  subspace  of  irreducible 
representations  of  S„t  spanned  by  the  outer  product  of  all 
antisymmetric  atomic  eigenstates, 

<J>(l;2;...;n:R)-«l>(l,;2';...;n':R)t 

N  N 

“=1  “=1  (16) 

where  i')  refers  to  the  Dirac  <5 function  in  the  totally 

antisymmetric  physical  subspace  (i:  Rtt)  of  the  atomic 
states  [cf.,  Eq.  (14)]. 

Although  some  of  the  non-totally-antisymmetric  repre¬ 
sentations  of  S„t  present  in  the  outer-product  basis  are  gen¬ 
erally  absent  in  the  basis  of  Eq.  (15),  a  certain  number  of 
these  non-Pauli  representations  are  always  present.  The  num¬ 
ber  and  type  of  unphysical  irreducible  representations  of 
spanned  by  the  reducible  spectral-product  representation  of 
Eq.  (15)  is  determined  by  the  total  number  of  electrons  n, 
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and  by  the  partitioning  of  these  into  the  particular  atomic 
values  >i\,  n2,-..,nN  employed.  Although,  strictly  speaking, 
these  unphysical  irreducible  representations  are  not  explicitly 
required  in  the  present  development,  they  can  be  enumerated 
in  specific  cases  employing  the  outer-product  reduction  rules 
in  the  usual  way.60  By  this  identification  it  is  made  clear  that 
the  basis  of  Eq.  (15)  can  span  unphysical  states  of  very  low 
energy  arising  from  “overfilled”  inner  shells  and  their  asso¬ 
ciated  continua,  in  which  the  physical  states  can  be 
embedded.32  37 

The  basis  of  Eq.  (15)  has  been  employed  previously  in  a 
number  of  related  connections,  including  long-range  and  ex¬ 
change  perturbation  theories25-38  and  so-  called  atoms- 
in-molec ules46  50  and  related  diatomics-in-molecules51-58 
methods.  In  contrast  to  the  present  development,  the  atoms- 
in-molecules  approaches  employ  explicit  prior  antisymmetri- 
zation  of  the  spectral-product  basis  of  Eq.  (15),  in  spite  of  the 
consequent  linear  dependence  in  the  limit  of  closure  this  ap¬ 
proach  entails,  and  the  perturbation-theoretical  methods  are 
generally  applicable  only  in  cases  of  weak  long-range  atomic 
interactions.  The  representation  of  Eq.  (15)  would  appear  to 
be  particularly  appropriate  for  molecules  and  other  atomic 
aggregates  even  in  the  absence  of  prior  aggregate  antisym¬ 
metry  in  that  it  is  strictly  orthonormal  and  contains  the  anti¬ 
symmetric  representation  only  once,  avoiding  thereby  prob¬ 
lems  of  linear  dependence.  It,  furthermore,  isolates  the  most 
singular  intra-atomic  Coulombic  interactions  in  the  indi¬ 
vidual  atomic  wave  functions  employed  as  a  basis,  gives  the 
exact  aggregate  energies  in  complete  dissociation  limits,  and 
results  in  a  Hamiltonian  matrix  of  particularly  simple  form 
(see  below).  Additionally,  the  basis  of  Eq.  (15)  provides  the 
possibility  of  developing  an  atomic-interaction-based  ap¬ 
proach  for  evaluating  an  aggregate  eigenspectrum  in  which 
the  exclusion  principle  among  aggregate  electrons  is  accom¬ 
modated  in  accordance  with  the  degree  of  interaction  be¬ 
tween  atoms,  rather  than  as  a  global  constraint  applied  uni¬ 
formly  and  uncritically  among  all  electrons  present,  which 
constraint  is  generally  not  required  in  the  case  of  highly 
separated  or  otherwise  weakly  interacting  atomic 

64 

constituents. 

IV.  SPECTRAL-PRODUCT  REPRESENTATION 
OF  SCHRODINGER  EIGENSTATES 

The  spectral-product  basis  defined  in  the  preceding  sec¬ 
tion  can  be  employed  to  construct  solutions  of  Eq.  (1)  fol¬ 
lowing  standard  variational  procedures.  Specifically,  employ¬ 
ing  an  expansion  of  eigenfunctions  in  the  spectral-product 
basis  of  Eq.  (15)  in  Eq.  (1)  and  projecting  with  each  of  the 
orthonormal  product  functions  in  the  basis  in  the  usual  way 
gives  the  familiar  matrix  Schrodinger  equation63 

H(R)-U*(R)  =  Uff(R)-E(R),  (17) 

where  E(R)  is  the  desired  diagonal  matrix  of  aggregate  ei¬ 
genvalues,  and  the  unitary  matrix 

UH(R)  =  (4»(l;2;...;n:R)|^(l;2;...;n:R))  (18) 

provides  a  representation  of  the  Schrodinger  eigenstates 
spanned  by  the  basis  of  Eq.  (15)  in  the  form 


vP(l;2;...;n:R)  =  <J>(l;2;...;n:R)  •  UH(R).  (19) 

The  semicolons  in  the  Schrodinger  eigenstates  *E(l;2;...;n:R) 
are  employed  to  emphasize  that  these  solutions  are  obtained 
from  Eqs.  (15)— (17),  in  which  the  orthonormal  aggregate 
product-function  basis  states  appearing  in  the  expression  of 
Eq.  (19)  for  the  wave  functions  are  not  individually  antisym¬ 
metric  under  all  electron-coordinate  permutations. 


A.  Form  of  the  Hamiltonian  matrix 

The  matrix  representative  of  the  aggregate  Hamiltonian 
operator  appearing  in  Eq.  (17)  takes  the  form45 

H(R)  =  («J>(  1;  2; . . .  ;n:  R)  ( 1, 2, . . . ,  n:  R)  |  <I>(  1;  2; . . . ;  n:  R)) 

N  f  N  -j 

=  2  H<“)+  2  V<“-«(R^)  ,  (20) 

a—  1  y  f3=  a+  1 

where  H(a)  and  V^a,^(Ra/3)  are  Hermitian  matrix  represen¬ 
tatives  of  the  operators  of  Eqs.  (3)  and  (4)  in  the  spectral- 
product  basis.  The  atomic  energy  matrix  H("'  is  a  constant 
diagonal  matrix  independent  of  atomic  position  vectors, 
whereas  the  interaction  matrix  V^a’^(Rap)  is  nondiagonal 
and  depends  explicitly  upon  the  vector  separation  Ray0  of  the 
two  indicated  atoms,  but  not  upon  the  individual  laboratory- 
frame  positions  Rfl  and  R^  of  the  two  atoms,  nor  upon  the 
position  vectors  of  the  other  atoms  in  the  aggregate.  The 
particularly  simple  form  of  Eq.  (20)  is  largely  a  consequence 
of  the  orthogonality  of  the  spectral-product  basis,  the  use  of 
atomic  eigenstates  in  the  representation,  and  the  atomic 
pairwise-additive  nature  of  the  interaction  terms  in  the 
Hamiltonian  operator  of  Eqs.  (2) -(4). 

The  specific  forms  of  the  atomic-energy  and  pairwise- 
atomic  interaction-energy  matrices  of  Eq.  (20)  depend  upon 
the  choice  of  ordering  rule  adopted  in  Eq.  (15).  Employing, 
for  example,  the  common  “odometer”  ordering  convention 
of  letting  later  indices  run  to  completion  prior  to  earlier  ones, 
the  last  atomic  matrix  H|A,)  and  the  pairwise  interaction- 
energy  matrix  \{N~1’N)(RN_UIV)  for  the  last  two  atoms  are 
seen  to  take  block  diagonal  forms,  with  matrices  constructed 
in  the  appropriate  atomic  and  atomic -pair  product  states,  re¬ 
spectively,  appearing  repeatedly  down  their  diagonals.  The 
atomic-energy  and  pairwise-atomic  interaction-energy  matri¬ 
ces  of  Eq.  (20)  for  an  arbitrary  pair  of  atoms  («,/3)  can  be 
obtained  in  this  convention  by  starting  from  the  block  diag¬ 
onal  versions  of  these  matrices  and  performing  appropriate 
row  and  column  interchanges  which  are  determined  by  the 
positions  of  a  and  /3  relative  to  the  last  two  indices  N—  1  and 
N,  respectively,  enumerating  the  atoms. 

Regardless  of  the  choice  of  ordering  convention  implied 
by  the  label  O  of  Eq.  (15),  the  atomic  matrices  can  always  be 
written  as  the  ordered  outer  products, 

H(“)  =  {ljj1)®l()2)«g)---®E^)®---®I<)Ar)}0,  (21) 

where 

Ka)  =  R  J I  '  a)(i:  R  J I  <Ha)(i:  R  J )  (22) 
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is  the  diagonal  matrix  of  physical  eigenenergies  of  the  atom 
a  [Eqs.  (9) — (11)],  arranged  by  convention  in  energy- 
increasing  order  down  the  diagonal.  The  interaction-energy 
matrices  are  correspondingly  written 

V(“’^(R  aP) 

=  {I< 1  >  ®  if  0  ■  ■  ■  ®  v<“’«(Ra/J)  ®  ■ ■  •  ■ ■  ®  if % ,  (23) 

where 
v(“'^(R^) 

=  <®(“-«(i;j:Rap)|V(“-«(l;j:Ra/J)|®(“-«(i;j:Ra/l)) 

(24) 

is  the  interaction-energy  matrix  in  the  ordered  atomic- 
product  basis  for  the  (a,/3)  pair,  given  by  the  expression 

<1>(  “•*>(!;  j :  Rff/s)  =  *>(i:  R  J  ®  j :  R*)}0 .  (25) 

In  Eqs.  (21)-(25),  if  is  the  unit  matrix  expressing  the 
orthonormality  of  the  antisymmetric  atomic  spectral  basis 
[Eq.  (9)],  an  outer-matrix-product  convention  is  adopted  for 
the  multiplication  of  matrices  enclosed  in  the  brackets  {•  •  •}, 
and  the  subscript  “ O ”  implies  ordering  of  the  terms 
in  the  energy  matrices  HUt),  V^“’^(Rajg),  and  vectors 
<|)(“,0)(i;j:R  )  jn  accordance  with  the  positions  of  the  indi¬ 
vidual  functions  in  the  row  vector  of  Eq.  (15).  The  pair- 
product  states  of  Eq.  (25)  are  written  for  notational  simplic¬ 
ity  as  functions  of  the  vector  separation  R([/3 ,  rather  than  of 
the  two  position  vectors  Ra  and  R^,  in  accordance  with 
primary  interest  here  in  their  use  in  evaluation  of 

(RffyS)- 

As  should  be  clear  from  the  discussion  of  Sec.  Ill,  the 
functions  <&(“'^)(i;j:  R^)  span  the  totally  antisymmetric 
representation  of  S„  +n^  once  and  only  once,  but  also  span 
other  irreducible  representations  of  the  a,  /3  atomic-pair  per¬ 
mutation  group,  and  so  these  functions  do  not  warrant  a  sub¬ 
script  “p”  designating  a  physical  subspace  of  the  generally 
reducible  product  of  Eq.  (25).  It  should  also  be  clear  that  the 
matrix  Ef  of  Eq.  (22)  in  the  atomic  spectral  basis  a  has  a 
dimension  smaller  than  that  of  the  matrix  v(“,/3>(Rajg)  of  Eq. 
(24)  in  the  a,  f3  atomic-product  spectral  basis,  although  the 
matrices  H(Q°  and  V(“’^(Raj8)  of  Eqs.  (21)  and  (23),  respec¬ 
tively,  have  the  same  dimension  as  the  entire  spectral- 
product  basis  of  Eq.  (15). 

The  vector  dependence  of  v^“’^(Rajg)  upon  Raj8  indi¬ 
cated  in  Eq.  (24)  arises  from  the  dependence  of  the  product 
basis  of  Eq.  (25)  on  both  the  magnitude  and  the  angular 
orientation  of  the  vector  Raj8  in  the  laboratory  frame.  Using 
the  transformation  properties  of  the  atomic  basis  functions 
under  rotations,65  the  interaction-energy  matrix  of  Eq.  (24)  is 
given  by  the  expression54 

v(“’«(R^)  =  R(“,/3)(Ra.j8)t-  y(a'p\Rap)  •  R^(RaP), 

(26) 

where  R(“^)(Raj8)  =  {D(“)(R„i8)®D(^(Rai8)}0  is  an  or¬ 
dered  outer  product  of  Wigner  rotation  matrices  for  each 
atom  effecting  the  transformation  from  a  reference  orienta¬ 
tion  in  which  the  two  z  axes  of  the  interacting  atoms  are 
co-aligned  to  one  in  which  they  take  their  physical  orienta¬ 


tions  in  the  laboratory  frame.65  The  irreducible  interaction- 
energy  matrix  ap)  of  Eq.  (26),  which  is  obtained 

from  Eq.  (24)  evaluated  in  the  co-aligned  reference  frame, 
depends  only  upon  the  scalar  separation  of  the  two  atoms 
Rap ,  whereas  the  notation  Ra/S(  =  4>ap ,  0ap$))  employed  in 
the  Wigner  functions  in  Eq.  (26)  refers  to  the  two  angles 
specifying  the  direction  of  the  vector  R,[/3  in  the  laboratory 
frame. 

Once  the  matrices  of  Eqs.  (22)  and  (24)  have  been  con¬ 
structed,  those  of  Eqs.  (21)  and  (23)  must  be  assembled  by 
performing  the  row  and  column  reorderings  implied  by  the 
symbol  “ O .”  Although  the  dimension  of  the  resulting  Hamil¬ 
tonian  matrix  of  Eq.  (20)  grows  exponentially  with  the  num¬ 
ber  of  atoms  considered  \^0(BN)  for  large  numbers  N  of 
identical  atoms  with  B  basis  states  each],  the  pairwise  nature 
of  the  interaction  terms  of  Eq.  (4)  ensures  that  nonzero  ma¬ 
trix  elements  are  obtained  only  between  aggregate  spectral- 
product  states  that  differ  by  no  more  than  two  atomic  eigen¬ 
states,  not  unlike  the  familiar  Slater’s  rules  for  two-electron 
atomic  matrix  elements  between  determinantal  wave  func¬ 
tions  comprised  of  orthonormal  orbitals.63  66  Accordingly, 
H(R)  is  a  sparse  matrix  which  contains  nonzero  elements  the 
number  of  which  in  any  row  or  column  grows  only  as  the 
square  of  the  number  of  atoms  N  and  of  the  number  of  basis 
functions  B  [«* 0(N2B 2)],  rather  than  exponentially,  with 
large  numbers  of  atoms  and  basis  functions.  The  latter  obser¬ 
vation  follows  from  the  number  of  spectral-product  functions 
that  can  differ  from  a  given  spectral  product  by  zero,  one,  or 
two  individual  atomic  labels,  and  ensures  that  H(R)  grows 
more  sparse  as  N  and  B  increase.  These  circumstances,  and 
the  fact  that  the  locations  of  the  nonzero  elements  in  any  row 
or  column  of  the  matrix  of  Eq.  (23)  are  known  a  priori, 
suggest  that  techniques  familiar  from  large-scale 
configuration-interaction  methods,63  including  particularly 
recursive  or  unitary  group  methods  which  avoid  construction 
of  the  entire  Hamiltonian  matrix  at  any  one  time,  can  be 
employed  in  obtaining  the  eigensolutions  of  Eq.  (17),  issues 
addressed  separately  below. 

B.  Evaluation  of  atomic-energy 
and  interaction-energy  matrices 

The  atomic-energy  matrices  of  Eq.  (22)  can  be  evaluated 
once  and  for  all  employing  conventional  electronic  structure 
theory  in  its  various  forms59-63  in  constructing  the  large  de¬ 
numerable  spectrum  of  eigenfunctions  generally  required  to 
satisfy  the  completeness  conditions  of  Eq.  (16).  Of  course, 
highly  accurate  experimental  values  are  also  generally  avail¬ 
able  for  a  significant  number  of  the  low-lying  atomic  ener¬ 
gies  of  most  atoms.67  Accordingly,  adopting  a  commonly 
employed  stratagem,  the  available  experimental  values,  pos¬ 
sibly  corrected  to  Coulombic  limits  by  approximate  removal 
of  relativistic  contributions,68’69  can  be  employed  in  place  of 
the  calculated  low-lying  values,  which  are  generally  less  ac¬ 
curate  than  the  former,  whereas  the  calculated  higher-lying 
values  can  be  employed  as  a  pseudospectrum  of  energies  in 
place  of  the  actual  Rydberg  and  continuum  atomic  states.70 
Such  a  calibration  procedure,46  which  is  easily  incorporated 
into  the  Hamiltonian  of  Eqs.  (20)-(26),  is  particularly  re- 


Downloaded  01  Nov  2007  to  198.14.54.253.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://jcp.aip.org/jcp/copyright.jsp 


9330 


J.  Chem.  Phys.,  Vol.  121,  No.  19,  15  November  2004 


Langhoff  et  al. 


quired  in  cases  where  the  calculated  low-lying  atomic  ener¬ 
gies  are  not  in  accord  with  the  energy  orderings  of  the  ex¬ 
perimental  spectra.  Employing  this  atomic  calibration,  the 
potential  energy  surfaces  obtained  from  the  development  are 
guaranteed  to  dissociate  into  accurate  atomic  separation  lim¬ 
its,  with  calculations  required  only  to  provide  (Coulombic) 
interaction  energies  among  the  atomic  constituents,  rather 
than  the  absolute  values  of  the  total  aggregate  energies.  Of 
course,  in  many  cases  such  an  atomic  calibration  procedure 
will  be  unnecessary,  and  the  calculated  atomic  energies  can 
be  employed  without  modification  directly  in  construction  of 
the  Hamiltonian  matrix. 

The  intergrals  required  in  the  irreducible  pairwise 
interaction-energy  matrices  of  Eq.  (26)  are  seen  to  be  two- 
center,  one-  and  two-electron  Coulombic  integrals,  which 
can  all  be  calculated  in  closed  analytical  forms  when 
Gaussian39  44  or  Slater-type71'72  orbitals  are  employed,  en¬ 
suring  that  high  precision  can  be  achieved  in  their  evalua¬ 
tion.  Specifically,  the  irreducible  interaction-energy  matrix 
appearing  in  Eq.  (26)  is  written  in  the  form 

y(a^\RaP) 


=  [  dra\  dvf{^a){ra)®^p\rp)}V{a’p) 
X(ra;r B-.Raf3) 

using  Eq.  (24),  where  [cf.,  Eq.  (4)] 

Zpe 2 


(27) 


V(“’w(r„;r  p:Rap)=- 


Zne 


R 


a/3 


K-Rsl  l^-RJ 


+ 


r„-r 


'p\ 


(28) 


is  the  two-center  Coulombic  interaction  operator  in  the  co¬ 
aligned  coordinate  system,  Ra=  R  ttk  and  R/3=/?^k,  with  k 
the  unit  vector  along  the  co-aligned  z  axes,  and 


y“)(0  =  (<“)(i:RJ|E  5(ra-r;)|<a)(i:RJ>  (29) 


is  the  so-called  one-electron  transition  density  matrix  for 
atom  a,73-75  with  a  similar  expression  employed  for  atom  f3. 

Evidently,  the  Hermitian  matrices  of  Eq.  (29),  properties 
of  which  have  been  discussed  previously  in  the  context  of 
nuclear  spectroscopy  and  nucleon  scattering  by  nuclei, 73-75 
determine  the  associated  irreducible  interaction-energy  ma¬ 
trices  of  Eqs.  (26)  and  (27),  and  so  constitute  the  computa¬ 
tional  invariants  required  in  the  development  of  Eqs.  (20)- 
(26).  The  matrices  of  Eq.  (29)  formally  involve  all  the 
eigenstates  of  an  individual  atom,  and  so  they  differ  from 
so-called  “reduced”  one-electron  density  matrices,  which  are 
scalar  functions  appropriate  for  individual  atomic 
eigenstates.62  Rather,  the  diagonal  elements  of  the  matrix 
y  “ 1  ( r  tt)  are  the  one-electron  charge-density  distributions  of 
the  spectral  eigenstates  of  atom  a,  whereas  the  off-diagonal 
elements  give  the  one-electron  transition  densities  connect¬ 
ing  these  eigenstates.  These  quantities  can  be  stored  in  the 
form  of  matrices  r(a)  in  the  combined  atomic-state  orbital- 
basis  (ra,ya)  representation,  and  employed  in  Eq.  (27)  in 


repeated  “on-the-fly”  calculations  of  the  interaction  energies 
required  in  construction  of  the  Hamiltonian  matrix  for  the 
atomic  configurations  R  considered.  Alternatively,  the  pair¬ 
wise  interaction-energy  matrices  of  Eq.  (27)  can  be  tabulated 
on  a  grid  of  interatomic  separations  Rap  for  the  interacting 
atomic  pairs  of  interest  and  retained  for  repeated  Hamil¬ 
tonian  assembly.  Accordingly,  in  this  approach,  large-scale 
configuration-interaction  calculations  of  the  atomic  spectral 
states  need  be  performed  only  once,  and  the  associated  one- 
electron  transition-density  or  interaction-energy  matrices  re¬ 
tained  for  repeated  use,  avoiding  explicit  high-level  molecu¬ 
lar  configuration-interaction  or  related  calculations  involving 
antisymmetric  aggregate  many-electron  basis  states  at  differ¬ 
ent  aggregate  configurations. 


C.  Eigenspectrum  of  the  spectral-product 
Hamiltonian  matrix 

Although  the  development  of  Eqs.  (1)— (29)  can  be  em¬ 
ployed  directly  in  computational  applications  in  selected 
cases,  a  number  of  factors  mitigate  against  adopting  such  an 
approach  more  generally.  Specifically,  the  spectrum  [Eq.  (1)] 
of  the  Hamiltonian  operator  of  Eqs.  (2)-(4)  contains  the 
aforementioned  unphysical  states,  including  continua  in 
which  the  physical  solutions  are  embedded  when  any  one 
atom  in  the  aggregate  contains  three  or  more  electrons.32-37 
Accordingly,  approximations  to  the  physical  states  of  Eqs. 
(l)-(4)  in  the  full  spectral-product  basis  may  generally  not 
be  easily  discernible  among  the  unphysical  quasicontinuum 
states  also  spanned  by  the  basis,  and  may,  in  fact,  be  com¬ 
putationally  inaccessible  in  the  presence  of  the  aforemen¬ 
tioned  low-energy  non-Pauli  states  present  which  are  associ¬ 
ated  with  overfilled  inner-shell  configurations  and  their 
corresponding  continua.  This  situation  is  complicated  by  the 
fact  that  precise  antisymmetry  will  only  be  obtained  from  the 
development  of  Eqs.  (1)— (29)  in  the  limit  of  spectral  closure, 
making  identification  of  the  desired  physical  states  poten¬ 
tially  subjective. 

The  task  of  classifying  the  states  obtained  from  the 
Hamiltonian  matrix  of  Eqs.  (20)-(29)  is  further  complicated 
by  their  degeneracies.  Specifically,  the  pairwise  interaction- 
energy  matrices  of  Eq.  (27)  will  be  diagonal  in  the  atomic 
spin  quantum  numbers  Sa,  Ms  and  S  p ,  Ms  ,  since  the 
Coulomb  interaction  operator  of  Eq.  (28)  is  independent  of 
spin  coordinates.  Consequently,  the  aggregate  Hamiltonian 
matrix  of  Eq.  (20)  will  generally  factor  into  noninteracting 
blocks  associated  with  particular  spin  configurations  of  the 
aggregate  atoms  ( Sa,Ms  ,a=  1,2,...,N),  and  the  eigensolu- 
tions  obtained  from  Eqs.  (17)-(29)  will  generally  be  frag¬ 
mented  into  degenerate  orthogonal  components  associated 
with  these  contributing  spin  configurations.  These  compo¬ 
nents  must  be  appropriately  combined  ex  post  facto  to  form 
states  of  the  correct  irreducible  symmetries. 

In  a  case  in  which  only  two  different  aggregate  spin 
configurations  are  present  in  the  spectral  product  basis 
[  «l>(  1;  2; . . . ;  n:  R)  =  <!>!  ( 1;  2; . . . ;  n:  R)  ©  <&2(  1;  2; . . . ;  N:  R)  ] ,  for 
example,  Eq.  (17)  factors  into  the  two  independent  equations 

H11(R)-U1(R)  =  U|(R)-E1(R),  (30a) 
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H22(R)  •  U2(R)  =  U2(R)  •  E2(R),  (30b) 

where  the  off-diagonal  matrix  elements  H12(R)  and  H21(R) 
vanish,  and  the  diagonal  energy  matrices  E^R),  E2(R)  and 
solution  matrices  U| (R),  U2(R)  are  obtained  from  the  indi¬ 
cated  separate  diagonalizations  of  Hn(R)  and  H22(R),  re¬ 
spectively.  The  eigenfunctions  of  Eq.  (19)  in  this  case, 

'E1(l;2;...;n:R)  =  <l>1(l;2;...;n:R)-U1(R)  (31a) 

tP'2(l;2;...;n:R)  =  <l>2(l;2;...;n:R)'U2(R),  (31b) 

may  include  a  solution  [''Pj(l;2;...;n:R)]t.  which  in  the  limit 
of  closure  has  the  same  eigenvalue  as  a  solution 
[^2(l;2;...;n:R)]/  if  both  spin  configurations  can  contribute 
to  the  total  wave  function  of  the  chemical  aggregate 
[{E1(R)}i.its“{E2(R)}//].  In  general,  these  two  components 
must  be  identified  and  appropriately  combined  to  constmct 
the  correct  antisymmetric  aggregate  eigenfunction  which  to¬ 
gether  they  comprise.  The  two  degenerate  spin  multiplets 
obtained  from  Eqs.  (31),  however,  need  not  necessarily  trans¬ 
form  irreducibly  under  electron  coordinate-label  permuta¬ 
tions,  and  so  there  may  be  no  way  in  which  to  readily  iden¬ 
tify  them  a  priori  as  potentially  degenerate  components  of  a 
totally  antisymmetric  physical  eigensolution  of  the 
Schrodinger  equation,  particularly  in  the  presence  of  other 
nearly  degenerate  unphysical  states.  Although  these  circum¬ 
stances  may  seem  unfamiliar,  a  similar  situation  would  arise 
in  treatments  of  atomic  structure  employing  unrestricted 
spin-orbital  product  configurations  in  constructing  the 
Hamiltonian  matrix,  rather  than  Slater  determinants.66  Com¬ 
binations  of  products  of  spatial  functions  and  spin 

1(-  •yn  HQ 

eigenstates,  ’  or  functions  of  spatial  coordinates  only, 
also  employed  in  constructing  eigensolutions  of  irreducible 
symmetry  similarly  correspond  to  solutions  expressed  in  the 
fragmented  forms  of  Eqs.  (30)  and  (31). 

Finally,  the  factoring  of  the  spectral-product  Hamiltonian 
matrix  into  individual  atomic  spin  blocks  can  also  lead  to  a 
failure  to  provide  the  correct  numbers  of  multiplet  states  in 
limits  of  large  interatomic  separations.  This  will  be  so  in  the 
common  case  of  dissociations  into  open-shell  atoms,  where 
the  failure  to  achieve  closure  in  the  spectral-product  basis  in 
the  absence  of  explicit  electron  exchange  terms  can  give  re¬ 
sults  which  include  only  the  lowest-lying  spin  multiplet  of  an 
allowable  number  of  degenerate  states  at  large  atomic  sepa¬ 
rations.  These  circumstances  can  lead  to  calculations  of 
highly  unphysical  potential  energy  surfaces  which  correlate 
with  incorrect  dissociation  limits.  Similar  remarks  apply  to 
cases  in  which  charge-transfer  terms  contribute  significantly 
to  eigensolutions  at  large  interatomic  separations,  where  the 
failure  to  achieve  closure  in  the  spectral-product  basis  and  in 
the  corresponding  development  of  Eqs.  (17)-(29)  could  in¬ 
troduce  spurious  variations  in  potential  energy  surfaces  with 
increasing  interatomic  separations. 


V.  CONVERGENCE  IN  THE  SPECTRAL-PRODUCT 
BASIS 

Group-theoretical  methods  have  been  used  previously  in 
understanding  the  consequences  of  the  change  in  electron- 
coordinate  permutation  symmetry  group  upon  formation  of 


atomic  aggregates,32  and  have  been  employed  above  to 
verify  the  presence  of  the  totally  antisymmetric  representa¬ 
tion  in  the  spectral-product  basis.  Such  methods,  however, 
are  largely  inapplicable  to  reduction  of  the  basis  in  the  group 
S„t  for  computational  purposes  in  light  of  the  nonsymmetri- 
cal  nature  of  the  representation  and  the  associated  absence  of 
all  electron  coordinate  permutations  in  explicit  form  in  Eq. 
(15).  In  the  present  section,  the  totally  antisymmetric  repre¬ 
sentation  of  S„t  spanned  by  the  basis  of  Eq.  (15)  is  isolated, 
the  additional  troublesome  issues  identified  in  the  immedi¬ 
ately  preceding  section  are  overcome,  and  convergence  to 
the  correct  Schrodinger  eigenstates  in  the  antisymmetric  sub¬ 
space  of  the  spectral-product  representation  is  demonstrated. 

Symmetry-adapted  eigensolutions  are  devised  in  the  fol¬ 
lowing  in  the  blocked  forms: 


1P(l;2;...;N:R) 


=  {ilr(p)(l;2;...;N:R),^(“)(l;2;...;N:R)}, 


E(R) 


E(f\R)  0  \ 

0  E<u)(R)) 


(32a) 

(32b) 


where  *P'(p)(l;2;...;n:R)  and  'E(")(l;2;...;n:R)  refer  to 
physical  “p”  and  unphysical  “u”  solutions  of  Eqs.  (17)- 
(19),  respectively,  and  E^R),  E^VR)  are  the  correspond¬ 
ing  diagonal  eigenvalue  matrices.  These  eigenspectra  are  ob¬ 
tained  formally  from  a  unitary  transformation  of  the 
aggregate  Hamiltonian  matrix  of  Eq.  (20), 


Hs(R)  =  Us(R)t-H(R)-U5(R) 

/H^(R)  0  \ 

\  o  h<u)(R)J’  (33) 

which  explicitly  isolates  the  physical  H$  (R)  and  unphysi¬ 
cal  (R)  blocks  of  the  spectral-product  Hamiltonian  ma¬ 
trix.  The  unitary  transformation  matrix  US(R)  obtained  from 
diagonalization  of  the  aggregate  metric  matrix  S(R)  con¬ 
structed  from  antisymmetrized  spectral-product  states  is 
shown  to  accomplish  the  blocking  of  Eqs.  (33),  and  the  cor¬ 
responding  eigenfunctions  and  eigenvalues  of  Eqs.  (32)  are 
shown  to  converge  in  the  limit  of  closure  to  results  obtained 
from  standard  variational  developments  employing  conven¬ 
tional  prior  antisymmetry  constraints. 

A.  Prior  antisymmetry 

The  development  leading  to  Eqs.  (32)  and  (33)  begins 
with  construction  of  an  antisymmetric  spectral-product  basis, 
familiar  from  the  early  development  of  Moffitt,46  obtained 
employing  the  prior  projection 


<I>(l,2,...,n:R)  =  HA(l,2,...,n)<I>(l;2;...;n:R),  (34) 


where  <t>(l;2;...;n:R)  can  be  taken  to  be  a  subspace  of  Eq. 
(15)  having  good  total  Ms  compatible  with  the  total  S  value 
of  interest  (Ms^S).  The  projector 


PA(l,2,...,n)  =  (»1!«2!---niv!)-1/2(u,!)-1/2E  (-1)^ 

p=  i 

(35) 
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is  the  familiar  n, -electron  antisymmetrizer  “wave-function 
normalized”  in  the  limit  R— >°°.63  In  Eq.  (35),  Sp  is  the  parity 
of  the  permutation  P p ,  the  sum  is  over  all  n , !  electron  co¬ 
ordinate  permutations,  and  the  correction  factor 
(«! \n2\'  ■  -nN\)~m  ensures  asymptotic  wave-function  nor¬ 
malization  in  the  presence  of  the  prior  antisymmetry  of  the 
individual  atomic  states  in  Eq.  (15).  The  use  of  commas  be¬ 
tween  electron  coordinates  on  the  left-hand  side  of  Eq.  (34) 
serves  to  indicate  that  the  individual  terms  in  the  basis  are  of 
good  permutation  symmetry  type,  as  opposed  to  those  of  the 
spectral  product  basis  of  Eq.  (15). 

Employing  Eqs.  (15)  and  (16),  the  basis  of  Eq.  (34)  is 
seen  to  satisfy  the  closure  relation 

0(l,2,...,n:R)-0(l',2',...,n':R)t 


n,\ 


nl\n2l---nN\  „=i  p 


n  4a)(i-i') 


(36) 


in  the  totally  antisymmetric  subspace  of  Eq.  (15).  Conse¬ 
quently,  the  antisymmetrized  spectral-product  or  Moffitt  ba¬ 
sis  is  formally  (2-fold  redundant  \Q  =  nt\l(nf.n2\---nN\)\ 
for  representations  of  antisymmetric  states.  ’  To  clarify 
the  origins  of  this  circumstance,  note  that  although  PA  of  Eq. 
(35)  nulls  the  unphysical  irreducible  representations  of  S„ 
contained  in  the  reducible  spectral-product  basis,  the  failure 
to  exclude  the  redundant  contributions  that  arise  in  antisym¬ 
metric  states  formed  from  unrestricted  orbital  products  in  the 
spectral-product  basis  gives  rise  to  the  factor  Q  in  Eq.  (36). 
Such  overcounting  is  avoided  within  each  individual  set  of 
atomic  states  ^'“'(iiRq.)  by  including  all  orderings  of  spin- 
orbital  configurations  but  excluding  explicit  permutations 
among  electron  coordinates  i,  or,  equivalently,  by  limiting 
the  contributing  spin-orbital  products  to  strictly  ordered  elec¬ 
tron  configurations  in  Slater  determinants.62  However,  since 
no  restrictions  are  placed  on  the  contributions  in  Eq.  (34) 
arising  from  products  of  spin  orbitals  from  two  or  more  dif¬ 
ferent  atomic  basis  sets,  and  all-electron  coordinate  permuta¬ 
tions  among  these  are  included  implicitly  in  the  limit  of 
closure,  Eq.  (34)  contains  redundant  terms  in  this 
limit.25,27-29,73-75  More  specifically,  the  factor  n,\  in  the  nu¬ 
merator  of  the  redundancy  factor  Q  is  the  total  number  of 
possible  arrangements  of  n,  electrons  among  n,  available  po¬ 
sitions,  whereas  the  factor  n  f.n2\  -  ■  -nN\  in  the  denominator 
of  Q  corrects  this  value  only  for  rearrangements  within  the 
individual  atoms,  which  do  not  arise  in  forming  the  linearly 
independent  atomic  spectral  states  by  using  a  single  assign¬ 
ment  of  electrons  to  orbitals  or  the  equivalent  restriction  to 
ordered  configurations  in  Slater  determinants.62 

The  overcompleteness  of  Eq.  (36)  is  of  no  consequence 
in  the  limit  R— since  the  interaction  terms  in  the  Hamil¬ 
tonian  matrix  vanish  identically  and  the  Hamiltonian  matrix 
correspondingly  has  no  off-diagonal  elements.  In  this  case, 
the  states  of  Eqs.  (15)  and  (34)  give  the  same  aggregate 
energies,  and  the  latter  otherwise  provides  results  equivalent 
to  use  of  the  former  atomic-product  states.  Moreover,  for 
large  values  of  R,  when  perturbation  theory  can  be  em¬ 
ployed,  the  factor  Q  can  be  taken  into  account  explicitly  in 
expressions  for  wave  functions  and  energies  involving  spec¬ 
tral  summations,  although  care  must  be  exercised  in  this 


connection.25-38  More  generally,  however,  for  arbitrary  finite 
values  of  R,  where  the  basis  of  Eq.  (34)  is  no  longer  ortho¬ 
normal  and  perturbation  theory  is  unsuitable,  the  redundancy 
in  the  limit  of  closure  gives  rise  to  linear  dependence  which 
must  be  overcome  in  computational  applications.  As  will  be 
demonstrated,  it  is  possible  to  not  only  employ  familiar  ca¬ 
nonical  orthogonalization  procedures  to  remove  the  linear 
dependence  in  Eq.  (34), 79  but,  correspondingly,  to  also  re¬ 
move  the  unwanted  non-Pauli  representations  of  the  symmet¬ 
ric  group  S„t  present  in  the  spectral-product  basis  of  Eq.  (15). 

B.  Metrically  defined  Hamiltonian  representation 

Following  standard  procedures,79  the  potentially  over¬ 
complete  Moffitt  basis  of  Eq.  (34)  is  employed  in  a  linear 
variational  solution  of  the  aggregate  Schrodinger  equation, 
giving  the  familiar  matrix  equation 

H(R)-Uh(R)  =  S(R)-Uh{R)-E(R)  (37) 

appropriate  in  a  nonorthonormal  representation.63  Here, 

N  (  N  -j 

H( R)=S  \Hia)( R)+  2  ^“’^(R)  (38) 

a=  1  y  f3=  a+  1 

is  the  metrically  defined  Hamiltonian  matrix  in  the  totally 
antisymmetric  basis,  where  the  individual  atomic  Hia)( R) 
and  atomic-interaction  V*“’^(R)  terms,  corresponding  to  the 
operators  of  Eqs.  (3)  and  (4),  respectively,  are  now  depen¬ 
dent  upon  the  positions  of  all  the  atoms  in  the  aggregate 
consequent  of  this  antisymmetry  [cf.,  Eq.  (20)], 

S(R)  =  (<!>(  1,2,..., n:R)|<I>(l,  2,..., n:R))  (39) 

is  the  overlap  or  metric  matrix  of  the  antisymmetric  basis, 
E(R)  is  the  aggregate  diagonal  energy  eigenvalue  matrix, 
and  UH(  R)  is  the  solution  matrix  that  contains  the  column 
eigenvector  representations  of  the  antisymmetric  Schro¬ 
dinger  states  in  the  nonorthogonal  basis  of  Eq.  (34), 

\P(l,2,  ...,n:R)  =  3»(l,2,  ...,n:R)  -UH(R).  (40) 

The  matrices  S(R)  and  UH(  R)  together  determine  the 
projection  of  the  basis  states  of  Eq.  (34)  on  the  solutions  of 
Eq.  (40)  in  the  usual  form, 

(0(l,2,...,n:R)|^(l,2,...,n:R))  =  S(R)  •  UH(R),  (41) 

where  UH(  R)  satisfies  the  familiar  modified  orthogonality 
condition 

Uh(R)'-S(R)-Uh(R)  =  I,  (42) 

and  I  is  the  unit  matrix  of  dimension  equal  to  that  of  the 
basis  of  Eq.  (34).  It  should  be  noted  that  different  fonts  are 
employed  in  Eqs.  (17)-(20)  and  in  Eqs.  (37)-(42)  to  distin¬ 
guish  between  matrices  and  other  quantities  represented  in 
the  spectral-product  and  Moffitt  basis  sets,  respectively. 

C.  Removal  of  linear  dependence  in  the  Moffitt  basis 

In  the  closure  limit  of  Eq.  (36),  the  Moffitt  basis  of  Eq. 
(34)  is  linearly  dependent,  the  metric  matrix  of  Eq.  (39)  cor¬ 
respondingly  becomes  singular  [det|S(R)|— >0],  and  Eq.  (37) 
is  generally  ill  conditioned.  The  linear  dependence  of  the 
basis  of  Eq.  (34)  in  the  limit  of  closure  is  avoided  by  isolat- 
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ing  a  linearly  independent  subspace  employing  the  unitary 
transformation  matrix  that  diagonalizes  S(R), 

Us(R)t  •  S(R)  •  US(R)  =  sd(R),  (43) 


where  the  diagonal  matrix  s£/(R)  contains  a  positive  block 
[s^(R)]  and  a  block  [s^u)(R)]  that  tends  to  zero  value  in 
the  closure  limit.  By  appropriate  orderings  of  the  eigenvalues 
and  vectors  of  S(R),  the  nonzero  block  associated  with  the 
column  eigenvectors  in  U$(R)  that  provide  the  physically 
significant  linearly  independent  subspace  of  Eq.  (34)  can  be 
placed  in  the  upper  left-hand  corner  of  s,/(R),  whereas  the 
block  that  tends  to  zero,  associated  with  the  linearly  depen¬ 
dent  or  “unwanted”  portion  of  the  basis  of  Eq.  (34),  can  be 
placed  in  the  lower  right-hand  corner,  giving 


srf(R) 


/  s^(R)  0  \ 

(Q*p  °) 

\  0  s<u)(R)r 

0  0  J 

(44) 


Here,  Ip  is  the  unit  matrix  in  the  linearly  independent  physi¬ 
cal  subspace,  0U  is  the  zero  matrix  in  the  unwanted  subspace, 
and  Q  is  the  redundancy  factor  of  Eq.  (36).  The  right-hand 
side  of  Eq.  (44)  follows  from  the  observation  that  the  eigen¬ 
values  of  the  metric  matrix  S(R)  are  either  Q  or  0  in  the  limit 
of  closure,  in  accordance  with  the  presence  of  this  factor  in 
Eq.  (36)  and  the  linear  dependence  of  the  unwanted  sub¬ 
space. 

Employing  the  matrix  U$(R)  of  Eq.  (43),  the  basis  of 
Eq.  (34)  is  transformed  into  a  new  basis  <55(1,2, ...,n:R) 
comprising  the  linearly  independent  subspace  <J>)f) 
(1,2,  ...,n:R)  in  which  the  physical  eigenstates  can  be  con¬ 
structed,  and  the  unwanted  subspace  <I>^(1,2, ...,n:R).  The 
two  subspaces  are  given  by  the  common  expression 

Os.(  l, 2, . . .  ,n:  R)  =  <5(  1,2, . . . , n:  R)  US(R)  sd(R)  “ 1/2 

=  {^)(l,2,...,n:R),^*)(l,2,...,n:R)}, 

(45) 

where 


<&s’)(l,2,...,n:R) 

=  {^(l,2,...,n:R)Us(R)}/,-s<f)(R)_1/2  (46) 

is  the  orthonormal  linearly  independent  subspace  and 
designates  the  physical  part  of  the  enclosed  row  vector.  The 
linearly  dependent  subspace 

1>f(l,2,...,n:R) 

=  {<5(1,2,. ..,n:R)-Us(R)}1(-s<u)(R)_1/2  (47) 

is  ill  defined  in  the  limit  sJ/'^R)— >0,. ,  but  is  not  required 
here.  The  closure  of  the  physical  basis  for  antisymmetric 
states, 

4>^)(l,2,...,n:R)-4>^’)(l',2',...,n':R)t 
=  {<5(l,2,...,n:R)-U5(R)}/,-s<f)(R)_1 

X  |<5(i',2', . . . ,  n' :  R)  •  US(R)}4— >  Q  ~  ‘{^(  l,  2, . . .  ,n:  R) 

X<5(i',2',...,n':R)}  — > j  U 

[  a=  1 


follows  from  Eqs.  (36)  and  (44)  and  the  unitarity  of  US(R), 
with  the  redundancy  factors  Q  that  appear  in  the  numerator 
and  denominator  of  Eq.  (48)  canceling  in  the  limit  of  com¬ 
pleteness.  The  second  subscript  “p”  on  the  right-hand  side  of 
Eq.  (48)  refers  to  restriction  of  the  ^-function  product  to  the 
linearly  independent  subspace. 

The  matrix  U$(R)  of  Eq.  (43)  can  be  employed  in  the 
usual  way  to  transform  the  potentially  ill-conditioned  Hamil¬ 
tonian  matrix  H( R)  of  Eqs.  (37)  and  (38)  to  a  form  that 
isolates  its  linearly  independent  block,  giving  the  generally 
well-conditioned  matrix  Schrodinger  equation 

H{f\  R)  •  l/^(R)  =  XfgX  R)  •  E(SPX  R),  (49) 

where  the  reduced-dimension  matrix 

H{sp\ R) 

=  <®^)(l,2,...,n:R)|ff(l,2,...,n:R)|^)(l,2,...,n:R)) 

=s(fXR)-m-{Vs{R)^H(R).\Js(R)}pp-S(fXR)-112 

^Q-m{Us(R)J,'-H(R)-Vs(R)}ppQ-m 

Nr  N  'j 

=  2  H{sa\R)+  E  V<“^(R)  (50) 

a=  1  y  /3—a+l 

provides  the  nonsingular  portion  of  the  metrically  defined 
Hamiltonian  matrix.  Here,  the  notation  {'"}pp  implies  that 
only  the  upper-left  {pp)  block  of  the  enclosed  matrix  is  to  be 
retained,  and  the  individual  atomic, 

Mff)(R) 

=  s<f)(R)-I/2-{U5(R)4.ff(“)(R).Us(R)};,p-s<f)(R)-1/2 

^e-1/2{Us(R)t.ff(“)(R).Us(R)}wg-1/2,  (51) 

and  atomic-interaction  terms, 

V^R)  =  s^(R)  - 1/2 •  {Us(R)t  •  V<a’«(R)  •  US(R)}W 
•s^(R)-1/2 

^Q- >/2{U,(R)t  ■  V<“-«(R)  •  US(R)}W(2- : m, 

(52) 

in  Eq.  (50)  refer  to  transformations  of  the  corresponding 
terms  in  Eq.  (38). 

The  matrix  of  physical  energies  R)  and  of  corre¬ 
sponding  column  eigenvectors  U^X R)  of  Eq.  (49)  are  of 
reduced  dimensions  relative  to  E( R)  and  UH(  R)  of  Eq.  (37), 
respectively.  Since  the  eigenstates  <5^(1, 2, ...,n:R)  of  the 
S(R)  matrix  of  Eq.  (41)  are  symmetry  adapted,  the  Hamil¬ 
tonian  matrix  of  Eq.  (50)  will  be  block  diagonal  in  the  total 
spin  quantum  numbers  and  in  other  appropriate  aggregate 
symmetry  labels,  in  accordance  with  Schur’s  lemmas.60  Of 
course,  the  limits  s^(R)— >QIp  and  s^")(R)-^0u  of  Eq.  (44) 
apply  only  when  closure  is  achieved,  with  the  partitioning  of 
sI/(R)  into  Sj  (R)  and  s{,u,(R)  somewhat  subjective,  and  the 
diagonal  matrix  (R)  appearing  in  Eqs.  (50)-(52)  required 
more  generally,  rather  than  its  limiting  value  QIp .  As  a  con¬ 
sequence,  selection  of  the  dimension  of  the  physical  block 
{’  ’ '  \ pp  °f  the  transformed  Hamiltonian  matrix,  which  will 
generally  have  nonzero  off-diagonal  blocks  connecting  the 
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“p"  and  “m”  subspaces  in  finite  representations,  is  corre- 
pondingly  subjective,  and  must  be  determined  in  the  course 
of  calculations  on  basis  of  the  accuracy  of  results  desired, 
with  convergence  achieved  when  the  off-diagonal  terms  are 
judged  to  be  negligible. 


D.  Isolation  of  the  antisymmetric 
spectral-product  subspace 

The  foregoing  standard  development  employed  to  re¬ 
move  linear  dependence  in  the  prior  antisymmetrized  or 
Moffitt  basis  of  Eq.  (34)  can  also  be  employed  to  separate 
the  physical  and  unphysical  subspaces  of  the  spectral- 
product  representation,  and  to  construct  corresponding  ener¬ 
gies  and  eigenfunctions  of  the  Hamiltonian  matrix  of  Eq. 
(17)  which  converge  to  the  symmetry-adapted  eigenstates 
obtained  from  the  development  of  Eqs.  (37)-(52).  Specifi¬ 
cally,  the  important  identity 

( <I>(  1;  2; . . . ;  n:  R)  |  _PA(  1, 2, . . . ,  n)  |  <&(  1;  2; . . . ;  n:  R) ) 

=  g-1/2S(R),  (53) 


which  follows  from  substitution  of  Eq.  (34)  into  Eq.  (39)  and 
is  valid  in  any  denumerable  basis,  indicates  that  the  matrix 
representative  of  the  antisymmetric  projector  of  Eq.  (35)  in 
the  spectral-product  basis  is  proportional  to  the  metric  matrix 
S(R)  constructed  in  the  antisymmetric  basis  0(l,2,...,n:R). 
Accordingly,  diagonalization  of  the  matrix  S(R)  in  Eq.  (43) 
is  equivalent  to  diagonalization  of  the  matrix  representative 
of  the  antisymmetric  projector  in  the  spectral-product  basis, 
and  to  the  construction  of  a  subspace  of  antisymmetric  eigen¬ 
states  in  this  representation. 

Applying  the  transformation  US(R)  to  Eq.  (53)  gives 

Us(  R) 1  •  ( <&(  1;  2; . . . ;  n:  R)  |  _PA  ( 1, 2, . . . ,  n)  1 0(  1;  2; . . . ;  n:  R) ) 


•us(R)=e 


,n(sAR> 


0 

s«(R) 


(54) 


where  Eq.  (43)  has  been  employed  on  the  right-hand  side. 
Defining  the  transformed  basis 


Equation  (57a)  demonstrates,  independent  of  group- 
theoretical  considerations,60  that  the  “p”  subspace  of  Eqs. 
(55)  and  (56)  contains  the  totally  antisymmetric  representa¬ 
tion  of  S„t  only  once  [sjf  ■'(Rj/Q— >Ip],  whereas  Eq.  (57b) 
demonstrates  that  the  “m”  subspace  consists  of  non-totally- 
antisymmetric  eigenstates  in  the  limit  of  closure  [s^“  (R)/<2 
— >0,.].  Note  that  the  factor  (1  /«,!)  employed  in  Eqs.  (57)  is 
appropriate  for  a  normalized  Young  operator,  rather  than  a 
normalized  wave  function  [cfi,  Eq.  (35)].  Finally,  from  Eqs. 
(16)  and  (56)  and  the  unitarity  of  US(R), 

^/“)(l;2;...;n:R)-a>^/H)(l';2';...;n':R)t 

II  ^(i— i')l  (58) 

“  1  >  plu 

expresses  the  closure  of  both  subspaces  of  the  eigenstates  of 
the  totally  antisymmetric  projector  in  the  spectral-product 
basis. 

It  is  seen  from  the  development  of  Eqs.  (53)— (58)  that 
both  “p”  and  “u”  subspaces  of  Eqs.  (55)  and  (56)  can  be 
constructed  employing  the  transformation  matrix  of  Eq.  (43) 
obtained  from  the  metric  matrix,  with  “p”  referring  to  the 
totally  antisymmetric  physical  subspace  of  the  spectral- 
product  basis  and  “u”  to  a  mathematically  well-defined  or¬ 
thonormal  unphysical  subspace  in  this  representation.  It  may 
be  concluded,  therefore,  that  the  transformation  matrix 
Uy(R)  of  Eq.  (43),  defined  to  eliminate  linear  dependence  in 
the  explicitly  antisymmetrized  Moffitt  basis  <I>(l,2,...,n:R)  of 
Eq.  (34),  correspondingly  separates  the  totally  antisymmetric 
and  non-Pauli  subspaces  of  the  spectral-product  basis 
<l>(l;2;...;n:R)  of  Eq.  (15).  Accordingly,  US(R)  can  be  em¬ 
ployed  to  partition  the  Hamiltonian  matrix  of  Eq.  (20)  into  its 
physical  and  unphysical  blocks  [Eq.  (33)],  and  to  construct 
physical  Schrodinger  eigenstates  in  the  spectral-product  ba¬ 
sis  [Eqs.  (32)]. 

Employing  the  transformation  matrix  UV(R)  of  Eq.  (43), 
the  physical  block  of  the  matrix  Schrodinger  equation  in  the 
spectral-product  basis  [Eq.  (17)]  takes  the  form  [cf.,  Eq. 
(49)] 


<&s(l;2;...;n:R)  =  <I>(l;2;...;n:R)  US(R) 

={^’)(l;2;...;n:R),  <^M)(l;2;...;n:R)}, 

(55) 

where 

^sP,U\  l;  2; . . . ;  n:  R)  =  {1>(  l;  2; . . . ;  n:  R)  •  Us(R)}p/u  (56) 

corresponds  to  the  physical  and  unphysical  subspaces,  it  fol¬ 
lows  that  Eq.  (54)  can  be  written  as  the  two  equations 


=  s(f](R)/Q^Ip,  (57a) 

j  «,! 

<0>'“)(l;2;...;n:R)|  —  2  (- l)^Pp|^“)(l;2;...;n:R)) 
nf  p=  i 

=  s',")(R)/e^0„.  (57b) 


H^(R)  •  U#>(R)  =  U^(R)  •  E<f  >(R),  (59) 


where 

Hf(R) 

=  (O^)(l;2;...;n:R)|ff(l,2,...,n:R)|O£0 
X(l;2;...;n:R))  =  {US(R)-H(R)-Us(R)}w 

N  (  N  } 

=  2  H<“;(R)+  2  V</-«(R)  (60) 

a=  1  y  (3=a+  1 

is  the  physical  block  of  the  Hamiltonian  matrix  of  Eq.  (20), 
with 

H^(R)={Us(R)t-H'“)-Us(R)}pp  (61) 

and 

Yj“’/3)(R)  =  {Us(R)1  •  V("^(Raj8)  •  Us(R)}pp ,  (62) 
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corresponding  to  transformations  of  the  individual  terms  of 
Eqs.  (21)  and  (23),  now  dependent  upon  the  positions  of  all 
the  atoms  in  the  aggregate  as  a  consequence  of  the  nonlocal 
effects  of  antisymmetry. 


E.  Equivalence  of  prior  and  post  antisymmetrization 

To  demonstrate  the  equivalence  of  results  obtained  from 
the  post  and  prior  antisymmetrization  procedures  described 
in  the  foregoing  developments,  note  from  Eq.  (34)  that  the 
states  <l>(/’)(l,2,...,n:R)  and  «I*^7)(l;2;...;n:R)  of  Eqs.  (46) 
and  (56)  are  identical  in  the  limit  of  closure, 

1,2, . . . ,  n:  R)  =  PA&/\l\ 2; . . . ;  n:  R)  -  Sj’)(R)  ~ 1/2 

nt  '■ 

(-l)^^)(l;2;...;n:R) 

p=  i 

— >  d>5P)  ( 1;  2; . . . ;  n:  R) ,  (63) 

where  the  left-hand  side  of  the  last  line  is  a  consequence  of 
Eq.  (44)  and  the  right-hand  side  of  the  antisymmetry  of  the 
basis  <f>5P)(l;2;...;n:R)  in  the  limit  of  closure,  the  n,\  per- 
mutations  of  the  un-normalized  antisymmetrizer  S  f=1 
(  — 1  )spP p  providing  n,\  identical  terms  which  cancel  the 
corresponding  factor  in  the  denominator  in  this  limit.  Alter¬ 
natively,  employing  Eqs.  (53)-(57),  the  spectral  overlap  of 
the  two  sets  of  basis  states  is  given  by  the  expression 

( ( 1;  2; . . . ;  n:  R)  |  <&^  ( 1, 2, . . . ,  n:  R) )  =  <2  1/2Sjp)(  R) 1/2 

-Ip,  (64) 

in  accordance  with  these  basis  states  becoming  term-by-term 
identical  functions  in  the  closure  limit  of  Eq.  (63). 

The  foregoing  results  can  also  be  expressed  in  terms  of 
the  representation  of  the  transformed  Moffitt  basis  [Eq.  (46)] 
in  the  transformed  spectral-product  basis  [Eq.  (56)], 

«I»s(l,2,...,n:R) 

=  <&s(l;2;...;n:R)  ■(<&s(l;2;...;ii:R)|<I>s(l,2,  ...,N:R)) 

- >  <&s(  i;2;  •  •  • ;  n:  R)  •  <2  -  1/2sd(  R)  ■  s^(  R) _  1/2 

-.*s(l;«...;i,:R).(1'  "J,  (65) 

employing  Eqs.  (53)-(57)  and  the  closure  limit  in  the  last 
line.  It  is  seen  that  Eq.  (65)  corresponds  to  Eq.  (63)  for  the 
physical  subspaces  of  the  two  basis  sets,  and  also  exhibits 
explicitly  the  linear  dependence  of  the  subspace  <J>!" * 
(1,2,  ...,n:R)  in  the  limit  of  closure. 

Similarly,  introducing  the  closure  relation  of  Eq.  (58)  for 
the  physical  states  on  either  side  of  the  Hamiltonian  operator 
in  Eq.  (50),  and  using  Eq.  (64),  gives 

H{sp)(  R) 

=  (^p)(l,2,...,n:R)|ff(l,2,...,n:R)|4»6,)(l,2,...,n:R)) 

=  Q~ms<d\ R) m- ■  s(f\R)  1/2g-1/2->H<p)(R), 

(66) 


employing  the  limiting  expression  s^lR)— >QIp  in  the  last 
line,  where  H'/YR)  is  given  by  Eq.  (60).  The  individual 
atomic, 

H(sa\  R)  =  Q  -  ms(f\R)  ~  m- H^(R)  •  s(dp)(  R)  “  mQ  “ 1/2 
-H<a)(R).  (67) 

and  pair-interaction, 

Vf’P 1  ( R)  =  Q  -  ms{f 1  ( R) +  1/2 •  V(sa/S)  ( R)  •  s(f 1  ( R) +  mQ  ~ 1/2 

->V^>(R),  (68) 

terms  of  Eq.  (50)  are  also  seen  to  be  equivalent  to  the  corre¬ 
sponding  terms  H^ff,(R)  and  R)  appearing  in  Eqs. 

(60) -(62)  in  the  limit  of  closure.  Equations  (66) -(68)  show 
that  the  Hamiltonian  matrices  of  Eqs.  (50)  and  (60),  and  the 
associated  individual  atomic-energy  and  atomic  pairwise 
interaction-energy  matrices  are  related  by  the  indicated 
renormalizations  in  finite  basis  sets,  whereas  they  are  identi¬ 
cal  in  the  closure  limit. 

Finally,  it  follows  from  the  foregoing  development  that 
the  eigenvalues  and  eigenvectors  obtained  from  Eq.  (59) 
converge  to  those  of  Eq.  (49)  [E^( R)— E^)(R),I7^)(R) 
— Uff^(R)],  alul'  in  view  of  the  equivalence  of  the  two  basis 
sets  [Eqs.  (63)-(65)],  that  the  eigenstates  obtained  from  the 
two  developments, 

■»P(p)(  1,2, . . . ,  n:  R)  =  1,2, . . . ,  n:  R)  •  U{^(R) ,  (69a) 

^(p)(l;2;...;n:R)  =  06,)(l;2;...;n:R)-U^)(R),  (69b) 

provide  identical  antisymmetric  solutions  of  the  Schrodinger 
equation  in  the  limit  of  closure. 

VI.  COMPUTATIONAL  IMPLEMENTATION 
AND  ILLUSTRATIVE  CALCULATIONS 

To  demonstrate  that  the  approach  described  in  the  fore¬ 
going  sections  is  computationally  viable,  illustrative  calcula¬ 
tions  of  the  lowest-lying  singlet  and  triplet  states  of  the  elec¬ 
tron  pair  bond  (H2)  are  given  as  an  example  of  the  nature  of 
the  metric  matrix  of  Eq.  (39)  and  its  eigenvalues,  and  of  the 
convergence  achieved  employing  standard  Gaussian  basis- 
set  methodology.39'40  A  computationally  efficient  implemen¬ 
tation  of  the  theory  is  described  which  employs  recursive 
construction80  81  of  invariant  subspaces  of  the  spectral  basis 
starting  from  appropriate  reference  states  in  calculations  of 
potential-energy  curves.  Construction  of  the  entire  metric 
matrix  is  avoided  in  this  way,  and  only  individual  rows  of  the 
spectral-product  Hamiltonian  matrix  need  be  evaluated  se¬ 
quentially,  avoiding  its  complete  construction  and  storage  at 
any  one  time. 

A.  Illustrative  calculations:  The  electron  pair  bond 

In  the  two-electron  case  (nt= 2),  the  spin  functions  are 
factored  out  and  attention  focused  on  spatial  functions  which 
are  symmetric  or  antisymmetric  under  electron  transposition. 
These  provide  a  basis  for  irreducible  representations  of  S2 
corresponding  to  single-  and  triplet-state  symmetries,  respec¬ 
tively,  with  the  singlet  states  playing  the  role  of  the  unphysi¬ 
cal  subspace  in  this  example.  That  is,  although  the  spatially 
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TABLE  I.  Spectral  energies  for  atomic  hydrogen. “ 


s  Basis*5 

p  Basis*5 

d  Basis'5 

/  Basis'5 

-0.499  991 

-0.124  998 

-0.055  256 

-0.026  707 

-0.124  994 

-0.053  439 

-0.011  199 

0.028  807 

-0.048  465 

0.025  060 

0.108  901 

0.187  597 

0.088  668 

0.274  825 

0.450  236 

0.617  389 

0.563  368 

0.986  192 

1.361  055 

1.748  440 

1.951  139 
5.660  465 
15.152  536 
38.982  050 
98.333  482 
249.224  098 
657.488  778 

2.880  038 

7.854  649 

21.314  359 

3.708  833 

9.719  458 
25.786  063 

4.791  877 

"Orbital  energies  (a.u.)  obtained  from  diagonalization  of  the  atomic  hydro¬ 
gen  Hamiltonian  employing  the  indicated  basis  sets  (Refs.  82-85). 
bBasis  sets  employed  are  the  most  diffuse  [  \2sH/)Hd 6/]  primitives  chosen 
from  a  set  of  14  even-tempered  Gaussian  primitives  for  each  angular  mo¬ 
mentum  symmetry,  which  set  consists  of  12  regularized  even-tempered 
primitives  (Ref.  82)  supplemented  with  two  additional  primitives  in  each 
symmetry  having  exponents  of  0.027  86  and  0.011  56. 


symmetric  states  are  physically  significant  for  H2 ,  they  nev¬ 
ertheless  represent  an  unwanted  symmetry  type  in  the  con¬ 
text  of  the  development  of  Sec.  V  relative  to  the  spatially 
antisymmetric  states,  and  so  can  be  used  to  illustrate  the 
permutation-symmetry  isolation  method  described  there.  In 
the  absence  of  other  irreducible  representations  for  n,  =  2, 
both  sets  of  solutions  can  also  be  constructed  by  direct  di¬ 
agonalization  of  the  spectral-product  Hamiltonian  matrix  for 
comparison  with  results  obtained  from  the  unitary  transfor¬ 
mation  of  Eqs.  (59)— (62)  for  isolating  the  subspace  represen¬ 
tations. 

The  spectral-product  basis  of  Eq.  (15)  in  this  case, 

®(l;2:/?)  =  {^“)(l:Ra)®^«(2:Ri8)}0,  (70) 

corresponds  formally  to  all  products  of  all  discrete  and  con¬ 
tinuum  hydrogenic  orbitals  for  the  two  atoms,  here  desig¬ 
nated  as  a  and  yS,  with  /?  =  |Ra—  Rg|  the  interatomic  separa¬ 
tion.  The  basis  of  Eq.  (70)  spans  both  spatially  symmetric 
and  antisymmetric  representations  of  the  group  S2  once  and 
only  once,  whereas  its  antisymmetric  form  [Eqs.  (34)  and 
(35)]  will  be  twofold  redundant,  as  is  seen  from  Eq.  (36).  To 
avoid  dealing  explicitly  with  continuum  hydrogenic  states, 
denumerable  representational  basis  sets  are  employed  in  the 
calculations  in  the  usual  way.39’40  Specifically,  approxima¬ 
tions  to  the  orbital  spectra  (/>'“'(  l:Ra)  and  2:Rg)  are 
devised  which  accurately  reproduce  the  lowest-lying  orbitals 
and  which  also  span  the  Rydberg  and  lower-energy  portion 
of  the  continuous  spectrum  in  the  form  of  pseudostates.70  By 
adopting  even-tempered  Gaussian  basis  sets  for  this 
purpose,82  standard  computational  methods  can  be  employed 
in  evaluating  the  matrix  elements  required  in  forming  the 
spectral-product  Hamiltonian  matrix  [Eqs.  (20)-(25)]  and  in 
performing  its  unitary  transformation  [Eqs.  (43)-(60)]. 39,40 

In  Table  I  is  shown  the  spectrum  of  atomic  energies  ob¬ 
tained  for  s,  p,  d,  and /orbitals  constructed  in  even-tempered 
Gaussian  basis  sets,82  with  orbital  exponents  chosen  so  that 
the  discrete  and  lower  continuum  states  are  spanned  by  the 
numbers  of  orbitals  shown  for  each  angular  momentum 


Langhoff  et  al. 


FIG.  1.  Eigenvalues  st  of  the  metric  matrix  of  Eqs.  (39)  and  (71)  for  H:, 
constructed  in  the  [12s8p]  basis  set  of  Table  I  as  functions  of  interatomic 
separation  R(a0)  (Refs.  84.  85),  values  ,s,~- 2  refer  to  approximately  anti¬ 
symmetric  eigenfunctions  of  electron  coordinates  constructed  in  the 
spectral-product  basis,  whereas  values  s,«*0  refer  to  approximately  symmet¬ 
ric  functions  of  electron  coordinates,  as  is  discussed  more  fully  in  the  text; 
apparent  discontinuities  in  the  slopes  of  the  plots  are  due  to  the  coarse  grid 
in  R  values  employed. 

value. 83-85  Although  systematic  studies  of  orbital  selection 
are  not  reported  here,  it  should  be  noted  that  considerable 
experience  has  been  gained  in  this  connection  through  pre¬ 
vious  studies  of  square-integrable  representations  of  all  the 
discrete  and  continuum  states  of  atoms  and  polyatomic  mol¬ 
ecules  constructed  on  a  common  basis.70  The  spectra  of 
Table  I  are  judged  on  basis  of  their  spatial  and  spectral  extent 
to  be  suitable  for  describing  the  charge  distortions  accompa¬ 
nying  chemical  bond  or  antibond  formation  in  H2  in  the  in¬ 
terval  R (1  -  5 )  a  Q ,  for  describing  wave  function  symmetry 
or  antisymmetry  in  the  absence  of  explicit  electron  exchange 
terms,  and  for  otherwise  approximating  spectral  closure  in 
this  interval. 

In  Fig.  1  are  shown  selected  eigenvalues  of  the  (1296 
X1296)  metric  matrix  S (R)  of  Eq.  (39)  for  H2  evaluated 
employing  only  the  s  and  p  basis  sets  of  Table  I.84’85  The 
normalized  two-electron  antisymmetrizer, 

Ra(1,2)  =  (1/2)-1/2(1-R12),  (71) 

is  employed  in  this  case,  with  the  redundancy  factor  Q 
=  (2  !/l !  1 ! )  =  2  indicating  that  the  totally  antisymmetric  rep¬ 
resentation  of  S2  is  spanned  twice  by  the  antisymmetric 
spectral-product  basis  in  the  limit  of  closure.  Only  the  fifty 
largest  (s,-®3 2)  and  the  fifty  smallest  (s,  =  0)  eigenvalues  of 
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S (R)  are  shown  in  Fig.  1  as  functions  of  the  interatomic 
separation  in  Hi.  The  eigenstates  of  S (R)  corresponding  to 
Sj*** 2  values  refer  to  linearly  independent  spatially  antisym¬ 
metric  states  constructed  in  the  antisymmetric  basis,  whereas 
those  corresponding  to  ,v ,  "=  0  values  refer  to  linearly  depen¬ 
dent  combinations  of  the  antisymmetric  basis,  in  accordance 
with  Eqs.  (43)-(48).  On  the  other  hand,  the  states  obtained  in 
the  spectral-product  basis  corresponding  to  the  sf*  2  of  val¬ 
ues  of  S(R)  refer  to  antisymmetric  states,  while  those  corre¬ 
sponding  to  5,^0  values  refer  to  symmetric  states,  in  accor¬ 
dance  with  Eqs.  (53)-(58).  When  the  two-electron 
symmetric  projector  is  employed  in  place  of  Eq.  (71),  results 
identical  to  the  foregoing  are  obtained,  but  with  the  states 
previously  corresponding  to  ,v  ,■  «=  2  and  ,v  ,■ 0  interchanging 
their  identities. 

A  significant  number  of  the  eigenvalues  of  S (R)  de¬ 
picted  in  Fig.  1  evidently  maintain  their  extreme  values  (s, 
=“0  or  2)  over  the  chemical  interaction  region  |  R 
s“(l-5)a0],  whereas  only  a  very  few  of  these  extend  into 
the  van  der  Waals  region  [/W“(5-  10)a0],  These  behaviors 
are  entirely  in  accord  with  the  spatial  characteristics  of  the 
spectral  states  of  Table  I  employed  in  constructing  S (R), 
which  states  have  relatively  small  amplitudes  at  distances 
*“(5-10)  «o  from  the  atomic  origins.  Accordingly,  the  basis 
of  Table  I  can  be  expected  to  give  converged  results  in  the 
chemical  interaction  region,  and  particularly  at  the  equilib¬ 
rium  interatomic  separation  (R=  1.40a0),  whereas,  alterna¬ 
tive  or  larger  basis  sets  are  required  to  achieve  closure  at 
larger  interatomic  separations. 

The  spectra  of  Table  I  are  employed  in  calculations  of 
energies  and  expectation  values  for  the  lowest-lying  'X*  and 

o  -j-  * 

X„  states  in  H2  at  the  equilibrium  interatomic  separation 
following  the  developments  of  Secs.  IV  and  V.84-85  These 
values  for  H2  can  be  obtained  directly  from  diagonalization 
of  the  full  spectral-product  Hamiltonian  of  Eq.  (20),  or  from 
the  unitary  transformation  of  Eq.  (33),  with  identical  values 
resulting  from  the  two  procedures  in  the  limit  of  closure.  The 
total  energies  and  binding  energies  and  the  expectation  val¬ 
ues  of  the  electron  transposition  operator  P 12  for  both  states 
shown  in  Table  II  evidently  converge  monotonically  in  both 
cases  to  known  values86  with  increase  in  the  angular  momen¬ 
tum  character  of  the  atomic  basis  states  employed.  Similarly, 
the  norms  of  the  familiar  singlet  and  triplet  Heitler-London 
functions87  represented  in  the  spectral-product  basis  also 
shown  in  Table  II  give  additional  indication  of  the  closure 
achieved  for  exchange  terms  in  the  wave  functions  in  this 
case.  It  is  found  in  the  finite  basis  sets  of  Table  I  that  signifi¬ 
cantly  more  than  one-half  of  the  transformed  Hamiltonian 
matrix  of  Eq.  (33)  is  required  in  order  to  reproduce  accu¬ 
rately  the  results  obtained  from  the  full  spectral-product 
Hamiltonian  matrix.  That  is,  the  symmetric  and  antisymmet¬ 
ric  subspaces  in  this  case  are  not  completely  isolated  by  the 
unitary  transformation  of  Eq.  (33)  into  approximately  one- 
half  of  the  dimension  of  the  original  spectral-product  basis, 
and  the  blocks  of  the  Hamiltonian  matrix  of  Eq.  (60)  identi¬ 
fied  as  {Uj(R)  •  H(R)  •  Us(R)}p/)/MU  must  include  a  sufficient 
number  of  small  but  nonzero  “off-diagonal”  contributions  in 
order  to  achieve  eigenvalues  in  accord  with  those  obtained 
from  the  full  Hamiltonian  matrix.  These  results  serve  to  em- 


TABLE  II.  Electron  Pair-Bond  Calculations.2 


Basis8 

Energy  (a.u.)c 

Binding  (eV)c 

(^12>C 

(|(p(HZ)|2)C 

H 

-1.0096 

l%g  state 

+0.2618 

+0.5255 

+0.8052 

[*/>] 

-1.0691 

+  1.8809 

+0.8409 

+0.9469 

[spd] 

-1.1140 

+3.1027 

+0.9525 

+0.9872 

[ spdf ] 

-1.1384 

+3.7667 

+0.9847 

+0.9967 

Exact  (Ref.  86) 

-1.1745 

+4.7478 

+  1.0000 

+  1.0000 

H 

-0.5586 

3£,t  state 

-12.0109 

-0.6226 

+0.2950 

[*/>] 

-0.6641 

-9.1400 

-0.6905 

+0.8078 

[spd] 

-0.7249 

-7.4856 

-0.9317 

+0.9538 

[spdf] 

-0.7524 

-6.7372 

-0.9801 

+0.9870 

Exact  (Ref.  86) 

-0.7842 

-5.8737 

-1.0000 

+  1.0000 

“Values  at  R—  1.40a0  obtained  from  diagonalization  of  the  Hamiltonian  ma¬ 
trix  of  Eqs.  (20)— (25)  for  H2 ,  or,  equivalently,  from  the  unitary  transfor¬ 
mation  of  Eq.  (33)  in  the  text  (Refs.  84,  85). 

bDenotes  the  portion  of  the  [  12i8p8d6/]  basis  set  indicated  in  Table  I 
employed  in  the  calculation. 

“Total  and  binding  energies  as  indicated;  (P  12)  refers  to  the  expectation 
value  of  the  electron  transposition  operator  P12 ;  |<p<Hi)|2  is  the  norm  of  the 
Heitler-London  function  as  represented  in  the  spectral-product  basis. 


phasize  that  the  development  of  Sec.  V  provides  a  formal 
proof  of  convergence  of  the  atomic  spectral  method,  rather 
than  an  optimal  computational  implementation  of  the  ap¬ 
proach,  and  that  additional  Gaussian  or  alternative  basis 
functions  and  other  computational  strategies  are  required  to 
improve  the  rate  and  degree  of  convergence  achieved. 


B.  Recursive  construction 
of  the  subspace  Hamiltonian 

The  metric  matrix  and  the  corresponding  explicit  unitary 
transformation  of  the  Hamiltonian  of  Sec.  V  are  avoided  em¬ 
ploying  a  recursive  method  to  generate  sequences  of 
(Krylov-Lanczos)  basis  states  which  span  the  two 
permutation-symmetry  subspaces  of  the  spectral-product  ba¬ 
sis  of  Eq.  (70).  The  development  employs  symmetric  or  an- 
tisymetric  starting  functions  in  its  implementation,  essen¬ 
tially  constructing  the  subspace  of  Eq.  (55) 

contained  in  Eq.  (70)  around  a  chosen  function  of  good  sym¬ 
metry,  and  avoiding  the  unwanted  subspace  (l;2:i?)  en¬ 
tirely.  In  this  way,  tridiagonal  Hamiltonian  matrices  of  much 
lower  dimensions  than  that  of  the  basis  of  Eq.  (70)  can  be 
constructed  which  converge  to  unitary  equivalent  representa¬ 
tions  of  the  totally  symmetric  or  antisymmetric  blocks  of 
H  (R)  without  explicit  evaluation  and  storage  of  the  entire 
metric  matrix  or  the  Hamiltonian  matrix  required  in  Sec.  V. 

Following  largely  standard  procedures,80'81  an  orthonor¬ 
mal  symmetric  or  antisymmetric  Krylov-Lanczos  subspace 
[Vj(l',2:R),j  =  l,2,...jmax]  of  the  spectral-product  basis  is 
written  in  the  form 

Vj(v,2:R)  =  <S>(l-,2:R)-Vj(R),  j=  l,2„..Jmax,  (72) 

where  the  column  vectors  Yj(R)  are  obtained  from  the  recur¬ 
sive  Lanczos  equations80'81 
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Pj(R)vj+1(R)  =  [H(R)- aj(R)I\  •  v/R ) 

-Pj-i(R)Vj-i(R),  j=  l,2,...j'max. 

(73) 

Here,  \o(R)  can  be  taken  as  a  null  vector  without  loss  of 
generality,  and  v \( R)  obtained  using  Eq.  (72)  from  a  suitably 
chosen  symmetric  or  antisymmetric  starting  function  (see  be¬ 
low).  The  coefficients  ctj(R)  and  f3j(R),  given  by  the  ex¬ 
pressions 


aj(R)  =  yj(R)l-K(R)-vj(R), 

(74) 

Pj(R)  =  vj+1(Ry-H{R)-Vj(R), 

(75) 

obtained  directly  from  Eq.  (73)  and  the  imposed  orthonor¬ 
mality  of  the  Vj(R),  provide  the  diagonal  and  off-diagonal 
elements  of  the  symmetric  or  antisymmetric  block  of  the 
Hamiltonian  matrix  in  tridiagonal  form.  A  sequence  of  ap¬ 
proximations  to  either  of  the  subspace  Hamiltonian  matrices 
is  obtained  in  this  way  with  increasing  order  ymax,  whereas 
convergence  is  achieved  in  the  limit  ymax  >'T-.7t) 

To  ensure  that  the  sequence  of  basis  states  of  Eqs.  (72)- 
(75)  generates  a  proper  subspace  of  the  spectral-product  ba¬ 
sis,  and  that  the  tridiagonal  Hamiltonian  matrix  of  Eqs.  (74) 
and  (75)  so  constructed  is  unitarily  equivalent  to  the  symmet¬ 
ric  or  antisymmetric  block  of  the  spectral-product  Hamil¬ 
tonian  matrix,  the  function  required  to  start  the  recurrence  of 
Eqs.  (73)-(75)  must  be  explicitly  symmetric  or  antisymmet¬ 
ric  and  have  nonzero  projection  on  the  entire  spectral- 
product  basis.  These  requirements  ensure  that  all  the  states  in 
the  sequence  of  Eq.  (72)  will  be  explicitly  symmetric  or  an¬ 
tisymmetric  in  the  closure  limit,  consequent  of  the 
permutation-symmetric  nature  of  the  Hamiltonian  in  Eq. 
(73). 

Once  the  Krylov -Lanczos  sequence  has  been  con¬ 
structed  employing  Eqs.  (72)— (75),  corresponding  Schro- 
dinger  eigenfunctions  are  obtained  in  the  form70 

J  max 

¥r(l;2:fl)=tfr(*)S  qj[Er(R)]v j(v,2:R) 

J  max 

=  NT(R)<t>(l-,2-.R)-lZ  qj[ET(R)]\j(R), 

7  =  1 

(76) 

where  Nr(R)  =  {'Zir™*qJ[Er(R)]2}-U2  is  a  normalization  fac¬ 
tor,  the  orthogonal  polynomials  q;{E)  satisfy  the  recurrence 
relation 

Pj(R)qj+i(E)  =  [E-aj(R)]qj(E) 

—  Pj-i(R)qj-i(E),  j=  i,2,...jmax, 

(77) 

and  the  starting  conditions  for  the  recurrence  are  qf]{E)  =  () 
and  q  ](E)  =  1 .  The  associated  energies  EV(R)  are  obtained 
from  the  secular  equation  qj  +1[£’r(/?)]=0,  with  the  qj(E) 
characteristic  polynomials  of  the  Hamiltonian  in  the  Krylov- 
Lanczos  basis.80'81  The  development  of  Eqs.  (73)-(77)  is 
identical  to  use  of  the  unitary  transformation  of  Eqs.  (53)- 
(62)  when  ymax  approaches  the  dimension  of  the  spectral- 


FIG.  2.  Total  energy  (a.u.)  curves  for  the  and  A  J  states  of  H2  as 
functions  of  atomic  separation  R(ct0 ):  heavy  solid  lines  refer  to  previously 
determined  accurate  values  (Ref.  86),  light  solid  lines  refer  to  Heitler- 
London  values  (Ref.  87),  dashed  lines  give  the  present  results  obtained  from 
the  recursion  procedure  indicated  in  the  text  employing  the  [s],  [jp],  and 
[sprf]  basis  states  of  Table  I  and  Heitler-London  test  functions  in  each  case 
(Refs.  83,  84). 


product  basis,  but  is  much  more  rapidly  convergent  to  spe¬ 
cific  portions  of  the  subspace  of  physical  solutions  than  is  the 
approach  described  in  Sec.  V.  This  is  because  appropriately 
chosen  starting  functions  v  \  (\\2:R)  can  act  as  points  of  spec¬ 
tral  concentration  about  which  optimal  Hilbert  space  descrip¬ 
tions  of  the  desired  physical  states  are  built  by  the  iterations 
of  Eqs.  (73)-(75). 70,80,81 

In  Fig.  2  are  shown  1 X ^  and  :i 2  „  potential  energy 
curves  in  H2  obtained  from  the  recursive  development  of 
Eqs.  (72) — (77)  and  the  basis  of  Table  I  employing  symmetric 
and  antisymmetric  Heitler-London  test  functions 
[u  j(l;2:/?)]  to  construct  the  starting  vector  v^/?)  in  each 
case.  These  reference  states  are  excellent  zeroth-order  ap¬ 
proximations  to  the  physical  states  desired,  with  the  leading 
tridiagonal  Hamiltonian  matrix  element  ( a- , )  the  Heitler- 
London  energy  in  each  case.  As  the  iteration  of  Eqs.  (73)- 
(75)  proceeds,  the  relatively  low-order  tridiagonal  matrices 
obtained  provide  rapidly  converging  approximations  to  the 
desired  subspace  energy  and  wave  function  corresponding  to 
the  chosen  starting  function,  avoiding  explicit  diagonaliza- 
tions  of  large-dimensioned  Hamiltonian  or  metrix  matrices. 
Evidently,  the  potential  curves  of  Fig.  2  obtained  converge 
rapidly  in  the  chemical  region  [f?  =  ( 1  -5)«o]  with  increase 
in  the  orbital  angular  momentum  values  employed,  the  [,v/;] 
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interaction  calculations  for  larger  R(3z\5a0)  values,  they  do 
not  reproduce  the  curvatures  in  the  configuration-interaction 
results  for  smaller  R  values.  This  circumstance  is  a  conse¬ 
quence  of  the  need  for  larger  or  alternative  basis  sets,  and 
possibly  higher  angular  momentum  terms  (see  Table  II),  re¬ 
quired  by  the  spectral  method  to  insure  convergence,  which 
terms  are  absent  in  the  results  of  Fig.  3.  That  is,  the  energy 
curve  in  the  \2a()-s  Rsi  I5a()  region  is  not  accurately  repre¬ 
sented  by  the  asymptotically  correct  Heitler-London  term 
and  the  [6s6/;>]  and  \6s6p6d]  spectral-product  basis  sets 
employed  in  Fig.  3,  emphasizing  the  well-known  sensitivity 
of  the  small  interaction  energies  in  the  long-range  region  to 
the  effects  of  electron  exchange.28'86 

VII.  DISCUSSION  AND  CONCLUDING  REMARKS 


FIG.  3.  Potential  energy  (a.u.)  curves  for  the  *2^  state  of  molecular  hydro¬ 
gen  in  the  long-range  interaction  region  obtained  as  indicated:  open 
squares — accurate  values  obtained  from  single  and  doubles  configuration- 
interaction  calculations  in  [6s6p]  and  f  6s6p6d~\  basis  sets  consisting  of  six 
regularized  even-tempered  primative  Gaussian  orbitals  of  each  angular  mo¬ 
mentum  symmetry  (Ref.  82),  solid  circles — spectral  theory  results  obtained 
from  the  recursive  development  of  Eqs.  (73)— (77),  calculated  employing  the 
indicated  basis  sets,  lines — fits  to  the  straight-line  portions  of  the  calculated 
values,  the  (1/R2)— >0  intercept  gives  C6  =  6.45  a.u.,  whereas  the  slope  of 
the  [6i6p6r/]  curve  in  this  limit  gives  Cs=  123  a.u.,  in  good  accord  with 
the  known  values  of  6.499  and  124.4  a.u.,  respectively  (Refs.  88.  89). 


limit  already  providing'-90%  of  the  chemical  bonding  en¬ 
ergy  at  the  equilibrium  interatomic  separation,  and  the  higher 
values  providing  results  in  increasingly  good  agreement  with 
the  correct  values.  The  Heitler-London  functions  for  H2  are 
well  represented  in  the  spectral-product  basis  [Eq.  (72)]  at 
the  equilibrium  internuclear  separation  [see  Table  II],  but 
must  be  included  explicitly  as  supplements  to  the  basis  at 
larger  separation  in  Fig.  2,  rather  than  employing  the  expan¬ 
sion  of  Eq.  (72).  In  this  latter  event,  one  additional  row  and 
column  is  added  to  the  matrix  H(/?)  in  Eq.  (73),  although  the 
development  of  Eqs.  (72)-(77)  is  otherwise  unchanged.84 

Finally,  although  the  present  implementation  of  the  de¬ 
velopment  was  not  devised  specifically  to  study  the  van  der 
Waals  region  of  bonding,  it  is  satisfying  that  the  correct  long- 
range  energy  limit  is  nevertheless  obtained  automatically 
from  the  recursive  projection  procedure.  In  Fig.  3  are  shown 
energy  curves  obtained  from  the  present  development  in  the 
long-range  region  (70=  12a0),  in  which  interval  the  singlet 
and  triplet  energies  are  indistinguishable,  in  comparison  with 
values  calculated  employing  standard  singles  and  doubles 
configuration-interaction  methods  in  the  indicated  basis 
sets.7,8  These  results  are  plotted  in  a  manner  [R6{E(R^>cc) 
—  E(R)}  vs  MR1]  that  allows  extraction  of  the  two  leading 
van  der  Waals  coefficients.  The  C6  coefficients  (»«6.45  a.u.) 
obtained  from  the  (l//?2)— >0  limits  of  both  sets  of  calcula¬ 
tions  in  the  [6s6/?]  basis  set  are  evidently  in  mutual  agree¬ 
ment  and  also  in  general  accord  with  the  known  correct 
value  (6.499  a.u.).88  Similarly,  the  C8  coefficients  («*123 
a.u.)  obtained  from  the  slopes  of  the  curves  in  Fig.  3  calcu¬ 
lated  in  the  \6s6p6d ]  basis  set  are  also  in  accord  with  the 
known  value  (124.4  a.u.).89  By  contrast,  although  the  energy 
curves  obtained  from  the  atomic  spectral  development  follow 
the  conventional  singles  and  doubles  configuration- 


The  atomic  spectral  method  described  here  adopts  a  dif¬ 
ferent  perspective  on  calculations  of  the  electronic  structures 
of  molecules  and  other  atomic  aggregates.  Enforcement  of 
wave  function  antisymmetry  in  an  atomic  spectral-product 
basis  is  deferred  until  after  construction  of  the  Hamiltonian 
matrix,  the  latter  consequently  taking  a  particularly  simple 
form  which  involves  only  unperturbed  atomic  energies  and 
pairwise-atomic  (Coulombic)  interactions.  The  required  in¬ 
teraction  energies  are  obtained  from  atomic  transition  density 
matrices  which  can  be  calculated  once  and  for  all,  providing 
thereby  a  quantitative  approach  to  the  quantum  theory  of 
matter  in  terms  of  spectral  attributes  of  its  atomic  constitu¬ 
ents.  Antisymmetry  requirements  are  subsequently  enforced 
employing  a  unitary  transformation  implemented  by  recur¬ 
sive  methods  devised  specifically  for  this  purpose.  These  and 
related  apspects  of  the  method  are  illustrated  with  calcula¬ 
tions  on  molecular  hydrogen  (H2),  which  demonstrate  con¬ 
vergence  to  results  obtained  from  conventional  prior  anti- 
symmetrization  methods.  Additional  applications  of  the 
formalism  presently  underway  should  help  to  demonstrate 
the  viability  of  the  development  in  cases  in  which  non-Pauli 
representations  are  present. 

The  theoretical  development  adopts  the  methodologies 
and  attributes  of  a  variety  of  conventional  methods  for  deter¬ 
mining  the  adiabatic  electronic  wave  functions  and  energies 
of  chemical  aggregates.  Use  of  the  atomic  spectral-product 
basis  dates  from  early  perturbation-theoretrical  treatments  of 
atomic  interactions23  and  related  later  efforts  to  incorporate 
the  effects  of  electron  exchange  in  this  context.26-38  In  con¬ 
trast  to  these  approaches,  the  present  development  is  nonper- 
turbative  in  nature  and  is  applicable  to  both  chemical  and 
physical  (noncovalent)  interatomic  interactions  on  a  common 
basis.  Although  explicit  electron  exchange  is  avoided  at  the 
outset,  the  totally  antisymmetric  subspace  of  the  atomic 
spectral-product  basis  is  isolated  and  the  unphysical  non- 
Pauli  representations  are  avoided.26-38  Convergence  is  en¬ 
sured  in  this  way  to  the  physically  significant  eigenstates 
contained  in  the  spectral-product  representation,  providing 
nonperturbative  methods  for  avoiding  the  non-Pauli  states 
that  accompany  the  change  in  symmetric  group  consequent 
of  aggregate  formation.32'37  Although  convergence  is  ensured 
in  the  limit  of  closure,  it  remains  to  illustrate  in  greater  detail 
the  computational  efficiency  of  the  method  described  for  re¬ 
cursive  isolation  of  the  physical  subspace  in  more  complex 
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cases  in  which  large  numbers  of  unphysical  discrete  and  con¬ 
tinuum  states  are  present.  Additionally,  although  no  explicit 
mention  is  made  here  of  perturbative  solutions  in  the  totally 
antisymmetric  subspace,  this  option  has  been  previously  con¬ 
sidered  and  remains  to  be  explored  in  complete  detail.45  Of 
particular  interest  in  this  connection  is  the  efficient  incorpo¬ 
ration  of  charge-transfer  and  other  explicitly  antisymmetric 
terms  in  the  form  of  chemically  significant  reference  func¬ 
tions  to  accelerate  convergence  of  the  development  to  anti¬ 
symmetric  states,  and  to  ensure  that  the  correct  numbers  of 
degenerate  multiplets  are  obtained  in  long-range  interaction 
limits. 

An  explicitly  antisymmetrized  atomic  spectral-product 
basis  has  been  employed  previously  in  so-called  atoms-in- 
molecules  approaches,46-50  dating  from  early  observations  of 
the  inadequacy  of  atomic-orbital  representations  in  describ¬ 
ing  atomic  multiplet  orderings  and  of  their  consequent  un¬ 
suitability  for  calculations  of  total  molecular  energies.46  The 
present  development  generalizes  the  notions  of  this  approach 
in  a  rigorous  manner  employing  the  atomic  charge  and  tran¬ 
sition  densities  of  the  noniteracting  atomic  constituents.  As  a 
consequence  of  Brillouin’s  theorem,90  such  one-electron 
properties  are  commonly  understood  to  be  less  sensitive  to 
the  effects  of  electron  correlation  than  are  total  molecular 
energies.  Once  the  atomic  charge-  and  transititon-density 
matrices  are  determined  employing  conventional  atomic- 
structure  methods,  they  can  be  employed  repeatedly  in  cal¬ 
culations  of  the  required  (Coulombic)  atomic  interactions. 
The  incorporation  of  experimental  atomic  energies  in  the 
Hamiltonian  matrix,  if  required,  is  straightforward,  and  the 
form  of  the  matrix  obtained  is  considerably  simpler  than  in 
the  original  atoms-in-molecules  approaches.46  50  Most  im¬ 
portant,  the  convergence  of  the  present  method  is  assured  in 
the  spectral  closure  limit,  and  difficulties  associated  with  po¬ 
tential  linear  dependence  are  avoided  entirely  in  the  ortho¬ 
normal  representation  employed.  In  light  of  these  attributes, 
largely  consequent  of  separation  of  Hamiltonian  matrix  con¬ 
struction  from  the  enforcement  of  wave  function  antisymme¬ 
try,  the  present  development  provides  an  atoms-in-molecules 
strategy  alternative  to  that  devised  many  years  ago  and  em¬ 
ployed  in  various  forms  in  the  interim.  Although  no  explicit 
mention  is  made  here  of  incorporating  diatomic  calculations 
or  other  related  information  in  aggregate  studies,51-58  this 
approach  has  been  described  previously  and  employed  in 
computational  applications  to  the  structures,  spectra,  and  dis¬ 
sociation  dynamics  of  atomic  clusters.45'91  The  important 
theorems  upon  which  these  studies  are  based  and  a  more 
detailed  accounting  of  this  variant  of  the  atomic  spectral  de¬ 
velopment  are  subjects  of  separate  reports. 

The  relatively  simple  forms  of  the  two-center  interaction 
terms  required  in  constructing  the  atomic  spectral-product 
Hamiltonian  matrix  allows  its  representation  in  a  wider  range 
of  orbital  functions  than  the  Gaussians  commonly  employed 
in  conventional  electronic  structure  calculations.  In  particu¬ 
lar,  Slater  forms  allow  for  efficient  evaluations  of  the  re¬ 
quired  interaction  integrals,71'72  potentially  providing  better 
representations  of  the  atomic  states  and  associated  one- 
electron  transition  densities  required  in  the  development  than 
those  obtained  in  Gaussian  representations.  The  ab  initio 


methods  currently  in  use  for  constructing  the  ground  and 
excited  electronic  states  of  atoms  and  molecules  can  be 
adopted  for  this  purpose  with  little  modification.7-12  Al¬ 
though  these  methods  as  currently  employed  can  provide 
highly  accurate  results  for  atomic  systems,  they  require  allo¬ 
cations  of  significant  computational  resources  in  order  to  per¬ 
form  the  repeated  electronic  structure  calculations  at  differ¬ 
ent  molecular  geometries  required  in  constructing  potential 
energy  surfaces  and  related  properties  of  molecules.  Addi¬ 
tionally,  they  generally  deal  with  total  electronic  energies, 
rather  than  with  atomic  binding  energies,  in  spite  of  the  well- 
known  difficulties  inherent  in  such  approaches,  the  constitu¬ 
ent  atoms  making  their  presence  apparent  only  through  the 
commonality  of  atomic  basis  sets  employed  in  such 
calculations.39-44  By  contrast,  the  present  approach  avoids 
repeated  high-level  molecular  calculations  of  total  aggregate 
energies  and  adopts  the  ab  initio  methods  only  in  accurate 
calculations  of  atomic  spectral  states  and  their  one-electron 
transition  densities.  Use  of  the  orthonormal  spectral-product 
basis  also  circumvents  the  well-known  basis-set  superposi¬ 
tion  difficulties  associated  with  conventional  approaches,92 
which  have  proved  particularly  troublesome  in  studies  of 
noncovalent  interactions.  Of  course,  systematic  methods 
must  be  devised  for  construction  and  tabulation  of  the  atomic 
information  required  for  implementation  of  the  present  ap¬ 
proach,  in  which  connection  it  should  be  noted  that  atomic 
spectra  are  of  considerable  interest  in  their  own  right,93  sepa¬ 
rate  from  their  use  in  the  context  of  calculations  of  aggregate 
electronic  structures. 

The  widely  employed  density-functional-theory  (DFT) 
approaches  provide  useful  structural  and  other  ground-state 
information  for  large  molecules  and  atomic  aggregates,13-17 
and  can  be  extended  to  studies  of  dynamical  response  func¬ 
tions  for  calculations  of  excited-state  properties.18  Although 
the  present  approach  is  based  on  largely  pedestrian  theoreti¬ 
cal  concepts,  and  entails  no  approximations  other  than  those 
of  a  computational  nature,  it  nevertheless  has  some  theoret¬ 
ical  commonality  with  the  DFT  methods  in  the  shared  at¬ 
tempt  to  circumvent  complexities  associated  with  antisym¬ 
metry  requirements,  although  by  different  means. 
Specifically,  conventional  antisymmetric  methods  are  em¬ 
ployed  in  calculations  of  the  required  atomic  spectral  infor¬ 
mation,  but  the  Coulombic  atomic  interaction  potentials  em¬ 
ployed  in  the  present  approach  to  aggregate  formation  are 
constucted  in  the  absence  of  inter-atomic  exchange  effects, 
and  are,  therefore,  largely  unphysical.  Relatedly,  attempts  to 
incorporate  explicitly  the  subtle  effects  of  electron  correla¬ 
tion  in  the  DFT  approaches  to  van  der  Waals  interactions 
apparently94  lead  naturally  to  the  one-electron  transition- 
density  matrices  of  particular  computational  focus  in  the 
present  development,  suggesting  that  study  of  possible  adop¬ 
tion  of  DFT  methods  in  the  present  context  may  be  benefi¬ 
cial. 

Finally,  although  it  is  premature  to  discuss  the  possible 
role  of  the  present  development  in  connection  with  semi- 
empirical  and  QM/MM  approaches  to  larger  systems,19-24 
there  may  ultimatly  be  some  relevance  in  this  connection  in 
view  of  the  atomic-interaction  perspective  adopted,  and  in 
light  of  the  growing  importance  of  such  methods  to  studies 
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of  proteins  and  other  biological  macromolecules.24  Particu¬ 
larly  needed  are  higher-level  QM  methods  to  treat  proton  and 
electron  transfer,  and  the  isomerizations  and  other  motions  of 
optical  chromophores,  in  the  presence  of  solvating  protein  or 
other  environments.  Additionally,  it  would  be  highly  desir¬ 
able  to  develop  methods  that  can  avoid  the  QM-MM  bound¬ 
ary  effects  experienced  in  some  applications  of  current  meth¬ 
ods,  among  the  other  challenges  such  studies  present. 

The  foregoing  issues  and  additional  matters  associated 
with  more  general  implementations  of  the  atomic  spectral 
approach  clearly  require  studies  well  beyond  the  scope  of  the 
present  paper.  Investigations  of  these  matters  and  the  results 
of  computational  applications  to  other  systems  currently  in 
progress  will  be  described  in  due  course. 
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